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Abstract. We discuss the dynamics of multiple scalar fields and the possibility of realistic 
inflation in the maximal gauged supergravity. In this paper, we address this problem in the 
framework of recently discovered 1-parameter deformation of SO(4,4) and SO(5,3) dyonic 
gaugings, for which the base point of the scalar manifold corresponds to an unstable de Sitter 
critical point. In the gauge-held frame where the embedding tensor takes the value in the 
sum of the 36 and 36 ’ representations of SL(8), we present a scheme that allows us to derive 
an analytic expression for the scalar potential. With the help of this formalism, we derive the 
full potential and gauge coupling functions in analytic forms for the SO(3) X SO(3)-invariant 
subsectors of SO(4,4) and SO(5, 3) gaugings, and argue that there exist no new critical points 
in addition to those discovered so far. For the SO(4,4) gauging, we also study the behavior of 
6-dimensional scalar fields in this sector near the Dall’Agata-Inverso de Sitter critical point 
at which the negative eigenvalue of the scalar mass square with the largest modulus goes to 
zero as the deformation parameter s approaches a critical value s c . We find that when the 
deformation parameter s is taken sufficiently close to the critical value, inflation lasts more 
than 60 e-folds even if the initial point of the inflaton allows an 0(0.1) deviation in Planck 
units from the Dall’Agata-Inverso critical point. It turns out that the spectral index n s of 
the curvature perturbation at the time of the 60 e-folding number is always about 0.96 and 
within the lcr range n s = 0.9639 ± 0.0047 obtained by Planck, irrespective of the value of the 
r] parameter at the critical saddle point. The tensor-scalar ratio predicted by this model is 
around 10 -3 and is close to the value in the Starobinsky model. 
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1 Introduction 

The recent observational data for CMB anisotropy [1] have definitely confirmed that the 
early universe underwent an accelerated expanding phase, dubbed inflation [2-4], Inflation 
is one of the most exciting topics in contemporary cosmology. Still, we do not have yet a 
satisfactory description of inflation from the viewpoint of fundamental theory, although a 
lot of endeavors have been dedicated since the proposal of the KKLT scenario [5, 6]. The 
most difficult issue to this problem is how to drive the cosmic acceleration while keeping all 
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moduli stabilized. The no-go theorem [7, 8] forbids the accelerating solution, provided we are 
discussing the problem within the framework of classical higher-dimensional supergravity (see 
also [9, 10] for the no-go theorem against stable de Sitter vacua in IIA theory). This theorem 
requires the necessity of introducing noncompact extra dimensions, higher-order corrections, 
non-perturbative effects and so on. 

Currently, not all of the lower-dimensional supergravities have been derived by dimen¬ 
sional reduction and dualities from the 11-dimensional supergravity. This fact motivates us 
to investigate inflationary models within 4-dimensional supergravity settings. Many people 
have focused attention to the N = 1 supergravity, since chiral fermions are naturally in¬ 
corporated in its theoretical framework. The N = 1 supergravity has been able to provide 
desired inflaton potentials by fine-tuning the (Hodge-)Kahler metric, the super potential, and 
possibly the D-term (see e.g, [11, 12] for a recent review). It is, however, difficult in general to 
embed these phenomenological models into string theory. Moreover, the N = 1 supergravity 
fails to have a strong predictive ability to explain the observational data, since the Kahler 
potential and superpotential are arbitrary functions. To the contrary, the potential for N > 2 
theories arises only from the gauging prescription so that the potential is highly restrictive, 
viz, we might be able to circumvent the landscape problem [13]. It deserves to comment that 
the implementation of this program is fairly nontrivial and requires an independent study, 
since N > 2 supergravities do not in general fit into the framework of the N = 1 super¬ 
gravity [14]. Extended gauged supergravities have attracted much attention recently also in 
the context of flux compactihcation and gauge/gravity correspondence. In particular, po¬ 
tentials arising from gaugings may have their origin in non-perturbative quantum effects in 
higher-dimensional theories. This expectation is supported by the fact that the potential is 
identified in some cases as (minus) the internal curvature of generalized geometry [15]. This 
prompts us to examine the possibility of inflation in extended supergravities. To this aim, the 
maximal N = 8 gauged supergravity certainly provides an optimal arena for exploring the 
gravity sector of string dynamics, since the IIA/IIB string theory and M-theory are formu¬ 
lated with maximal supercharges corresponding to IV = 8. Over the past 30 years, the N = 8 
gauged supergravity [16] has been intensively studied from various points of view [17-34], 1 A 
standard lore that had emerged from these analyses is that the maximal gauged supergravity 
is unique. This stereotyped idea was dispelled by the recent discovery of a new 1-parameter 
family of deformation of SO(8) gauged supergravity [36], based upon the embedding tensor 
formulation [37, 38]. Although the original SO(8) gauged supergravity has been obtained via 
gaugings and also by the S' reduction of the 11-dimensional supergravity [39], the uplifting 
of the deformed theory into the 11-dimensional framework is yet unknown (see [40] for a 
work along this direction). A distinctive feature of the new deformation is that it introduces 
“magnetic” vector potentials, in addition to the standard electric ones. The introduction of 
magnetic potentials manifests itself in a variety of ways. An interesting consequence of this 
dyonic gaugings is that they possess a completely distinct structure of extrema, compared 
to the standard electric ones. Combined with the “go to the origin” approach developed in 
refs. [41, 42], we can now obtain dozen of new critical points, new supersymmetry breaking 
patterns, the analytic mass spectrum for seventy scalars and a black hole solution [41, 43-53]. 
In our previous paper [43], we revealed that the special SL(8)-type critical points 2 can be 

1 See ref. [35] for a recent attempt to embed the standard model fermions in IV = 8 supergravity. 

2 We refer to a gauging with its gauge group contained in the standard SL(8) maximal subgroup of £ 7 ( 7 ) as a 
SL(8)-type gauging. If a critical point of the potential in a SL(8)-type gauging can be further transformed into 
the origin corresponding to the coset SU(8) in the scalar manifold U 7 ( 7 )/SU(8) by a duality transformation 
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completely classified and their mass spectrum can be evaluated group-theoretically. It turned 
out that many of SL( 8 )-type gaugings allow a deformation parameter as in the SO( 8 ) theory. 
An intriguing facet of special SL( 8 )-type (anti-)de Sitter extrema is that their mass spectrum 
is only dependent on the residual gauge symmetry, insensitive to the original gaugings and the 
deformation parameter. Further, special de Sitter critical points appear only for SO(4,4) and 
SO(5, 3) gaugings and are unstable in the directions that are singlet with respect to the resid¬ 
ual symmetry. The parameter r] = V"/V characterizing the instability growth rate is 77 = —2, 
which implies that the instability grows in a time scale comparable to the cosmic expansion 
time in the cosmological context. Remarkably, the adjustable deformation parameter can be 
absorbed in the overall cosmological constant and does not affect the r/ parameter. This is 
a typical r^-problem by which many of supergravity theories are plagued. Accordingly, a de 
Sitter critical point around which a realistic inflation can be derived must be of the general 
type sitting off the origin of the scalar manifold in any SL( 8 )-type gauge frame, if it exists. 

In ref. [46], Dall’Agata and Inverso found a new de Sitter critical point of general type in 
the SO(4,4) gauging (we shall refer to this extremum as the DI critical point) which evades 
the 77 -problem (see ref. [55] for the 11-dimensional uplift). The modulus of the 77 parameter 
at this critical point can be made as small as we wish by choosing the deformation parameter 
close to a critical value. In the present paper, we will make deeper studies of the scalar 
held dynamics in this gauged supergravity by extending the analysis of [46]. We present a 
simple method to compute the scalar potential away from the origin, based on an inventive 
diagonalization of the coset representative . Applying this prescription to the SO(3) x SO(3)- 
invariant sector of SO(4,4) and SO(5, 3) gaugings, we examine whether there exist additional 
critical points in these six-dimensional sectors or not. As far as the authors know, the present 
study provides the highest-dimensional analytic search of critical points in N = 8 gauged 
supergravity with noncompact gaugings. We further deeply investigate whether a realistic 
inflation is realized around the DI saddle point by solving the equations of motion in the 
six-dimensional SO(3) X SO(3)-invariant subspace of the 70-dimensional scalar manifold in 
the SO(4,4) gauging. We also take into account the flux contributions and examine if the 
chromo-natural inflation [56, 57] and the anisotropic inflation [59, 60] occur. We find that 
by choosing the deformation parameter close to the critical value, the universe can inflate 
for a sufficiently long time even if the initial offset of the inflaton from the DI saddle point 
is of order unity in the Planck units. We also find that the spectral index n s of the scalar 
curvature perturbation at the present horizon scale is quite insensitive to the initial condition 
and the existence of flux and consistent with the observations by Planck [1], We will also see 
that the tensor-scalar ratio predicted by our model is of order ICR 3 and about one third of 
the value predicted in the Starobinsky model [2]. 

Because the potential of this gauged maximal supergravity has no lower bounds, we 
cannot discuss the reheating phase and subsequent big-bang stage. In order to investigate 
these post inflation stages, we have to find a modification of the theory in such a way that 
the potential is non-negative and have to couple the system to a matter sector with a lower 
local supersymmetry. In the present paper we do not attempt to make such modifications, 
and examine the scalar dynamics by assuming that the inflationary stage is well described by 
the maximal gauged supergravity. 

The present paper is organized as follows. In the next section, after fixing our notations, 
we develop a general formulation that can be used to derive an explicit expression for the 
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potential in the SL( 8 )-type gauge frame. In sections 3 and 4, we apply this formulation to 
the SO(3) x SO(3)-invariant sector of the SO(4,4) and SO(5,3) gaugings, respectively. We 
derive explicit expressions for the kinetic term and potential for the scalar fields as well as the 
gauge coupling functions, and perform the search of new critical points in this sector. Next, 
we do cosmology in section 5. We first derive the effective action describing the cosmological 
evolution of the scalar-gravity-gauge system obtained in section 3. We clarify the definitions 
and meanings of slow-roll parameters in our multi-component system. In particular, we point 
out the existence of special attractor slow roll trajectories that play a crucial role in the 
inflation dynamics in our system. After analyzing the ?/-problem in our system in detail, 
we explore the possibility of a realistic inflation by solving the multi-dimensional dynamical 
evolution equations numerically. Section 6 summarizes the results of the present work. 


2 SL(8) gauge-field frame 


We begin this section with a brief review of N = 8 supergravity results in the existing 
literatures [38, 41, 43], based on the embedding tensor formalism. Our notations basically 
follow refs. [38, 41, 43]. 

The embedding tensor enables us to discuss the embedding of a gauge group G into 
the duality group E 7 ^ in a group-theoretical fashion [37, 38]. The embedding tensor ®M a 
specifies the gauge generators Xm hi terms of the E 7 ^ generators t a as Xm = ®M a t a - Here 
= 1,..., 56 and a, (3, ... = 1, . 133 are the fundamental and adjoint representation 
indices of E 7 ^ 7 y Using the gauge generators Xm, the gauge coupling of scalar fields can be 
introduced in the coset representation by the replacement d^ —y d^ — gA^ M X m, where g is a 
gauge coupling constant and A^ 1 is a gauge potential. 

Consistency of a gauging amounts to requiring that the embedding tensor obeys a set 
of linear and quadratic constraints. The linear relation renders ®M a to belong to the 912 
representation of £ 7 ( 7 ), which branches into 912 — > 36 + 36 7 + 420 + 420 7 under its maximal 
subgroup SL( 8 ). In this section, we restrict our consideration to SL( 8 )-type gaugings and 
work in the SL( 8 ) gauge-field frame. Then, only the 36 and 36 7 representations are turned 
on, and the gauge generators Xm = T {X ab X ab ) take the form, 


X ab = 


_ ( S[a e 9 b \[ c S d ]fl 0 


0 -fi[a C 9b][efif] d ^ 


x ab = 


-V a ^ 1[< V ] 0 


e [a a b][c 5 f] d] 


( 2 . 1 ) 


Here, 6 ab and £ ab are 8 x 8 symmetric matrices, and a,b,... = 1,...,8 are SL( 8 ) indices. The 
quadratic constraint requires the following relation [41] 


e • e = ci 8 , 


( 2 . 2 ) 


where c G K is a constant and I n is the n-dimensional unit matrix. The electric gaugings 
£ = 0 trivially satisfy the quadratic constraints. 

In maximal supergravity, scalar fields are described by the symmetric coset space -EV( 7 )/SU( 8 ). 
In gauged versions, these scalars acquire a potential to keep supersymmetries. In the embed¬ 
ding tensor formulation, the potential is also expressed in an £ 7 ( 7 ) covariant way and reads 

v=^- (X mn r X pq s M mp M n Qm rs + 7Xmn Q X pq n M mp ) , (2.3) 
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where Mmn = (£ ■ t L)mn is the 56 x 56 scalar matrix with its inverse (M -1 ) MAr = M mn = 
Q MP £l N Q Mpq. Here, & IN = Qmn = (*C 2 <S> l 28 )jwiv is an Sp(56, M)-invariant metric, and 
L = L((j)) is a coset representative in the Sp(56,M) representation described below in more 
detail. The kinetic term of the scalar fields can be expressed as e _1 £ = — p ijkiV^ kl , 
where V lu jki is the 77 -selfdual part of the T-tensor [see eq. (2.22) below for definition]. In 
the absence of the gauge held, this kinetic term can be regarded as the Lagrangian for the 
harmonic map from the spacetime to the target scalar manifold E 7 ^/SXJ(8) with the metric 

d4 = y/m ® r ijkl = -ifr (dMdM- 1 ). (2.4) 

We call this metric a kinetic-term metric or a scalar manifold metric in the present paper. 

In appearance, the potential (2.3) is a highly nonlinear function of 70 scalars, which is the 
main obstruction when one attempts to find critical points. A proposal put forward in [41, 42] 
is to move a critical point to the origin of the scalar manifold by an £ 7 ( 7 ) transformation acting 
naturally on the scalar manifold. This procedure transforms the embedding tensor as well. It 
follows that the critical point condition d p V = 0 gives rise to another quadratic condition on 
the embedding tensor [41]. Here p = 1,..., 70 runs over only the noncompact directions of the 
coset space £ 7 ( 7 )/SU( 8 ). This additional relation must be solved together with the quadratic 
constraints of the embedding tensor, and these algebraic equations are indeed amendable 
to analytic study. In our previous paper [43], we demonstrated that all special SL( 8 )-type 
critical points, which can be mapped to the origin by SL( 8 ) transformations, are exhausted 
by known ones." It turned out that only the SO(4,4) and SO(5,3) gaugings admit de Sitter 
critical points of the special type. Both of these critical points possess the tachyonic directions, 
irrespective of the deformation parameter, providing candidates of an inflaton. Unfortunately, 
the curvature of these critical points turned out to be too large to allow a sufficient duration 
of inflation. In the present paper, we explore their “off-center” behavior, which was outside 
the scope of our previous analysis [43]. 


2.1 Coset representative 

To find critical points away from the origin, we need to compute explicitly the exponential 
map of the coset representative. In this section, we present an efficient formulation that 
facilitates to perform this program for SL( 8 )-type gaugings. 

The 70 scalar fields take values in the 70 r representation of SU( 8 ) and are expressed as 

tfijkl 4*[ijkl\ A £ijklmnpq ( P ^ ■ (2.5) 

Here and in what follows, the SU( 8 ) indices i,j,... = 1,..., 8 are raised and lowered by complex 
conjugation, i.e., (fijki = 4> l ^ kl ■ In the Usp(56) representation [61], we can adopt the unitary 
gauge in which the £ 7 ( 7 ) coset representative 5? takes the form, 


5? = exp 


0 <A 

0) 


( cosh(y#t) ^)jh( v /^) \ 

^tjh(vW) cosh( v / 0 t ~4>)) 


E Usp(56), 


( 2 . 6 ) 


3 In ref. [43], it has been argued that the mass spectrum of the scalars is dependent only on the residual 
gauge symmetry, rather than the original gaugings when the cosmological constant is nonvanishing. We have 
also confirmed that this property is shared also by other fields with spin. 
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where we have defined 


jh(z) = z 1 sinh(z). 


(2.7) 


In the double-index notation, (ftijki = <t>[ij][ki] should be viewed as a 28 x 28 symmetric matrix. 
Noticing (f> jh(-y/ <j>^(j>) = jh(-y/the following properties hold 






( 2 . 8 ) 


Our remaining task is to find how to implement the exponential map for the coset 
representative (2.6). This can be tackled most easily as follows. Since cj) can be viewed as a 
28 x 28 complex symmetric matrix, the Autonne-Takagi factorization [62] enables us to find 
matrices A and U such that 


4> = UA t U , U € U(28). (2.9) 

Here A = diag(Ai,..., A 28 ) is the diagonal matrix with real entities. Equation (2.9) is tanta¬ 
mount to the following eigenvalue problem, 

<j)u = Xu , Uij = ( Uj)i . (2.10) 

If this system is solved, the square-root matrix can be easily computed to give yj 0</>t = UhXJ\ 
thereby the coset representative (2.6) now takes the simple form 

_ / U cosh A£/t U sinh A T f/ A . . 

■ = y U sinh AU^ t7coshA T t/y ' 

The computation of the scalar matrix Mmn asks for the change of the basis from the Usp(56) 
basis to the Sp(56,M) basis. This can be done through the Majorana-Weyl SO(8) gamma 
matrices [61], 

•sV 1 = Sj£ p rp°(S' 1 )«-, s ^ jA (A.“‘ jjg* J = (s*)- 1 . (2.12) 

Here Lm N is the coset representative in the Sp(56,R) representation and the underlined 
indices refer to the SU(8). In the Sp(56,M) basis, we have 

Y^( r ^ r )ab a/ = 2 + iU abcd , S ab = S( ab ), U abcd = U [ abcd ] = -kU abcd . (2.13) 

Each of real tensors (S ab , U abcd ) has 35 components. Since there is no need to distinguish the 
upper and lower indices of gamma matrices, we can identify from the scratch, 

—T^>r —>• i f>, S —>• <S> I 28 > (2-14) 

16 

where 

*=^ 0 -1) ■ (2i5) 

With this identification, the matrix Mmn reads 

Mmn = (S^ 2 S)mn = ( ) , (2-16) 
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where 


Mabcd =\[(U + U)e 2A T (U + U)-{U- U)e~ 2A T (U - U)} abcd , (2.17a) 

M ab cd = \[{U- U)e~ 2A T (U + U)-{U + U)e 2A T (U - U)\ ab cd , (2.17b) 

M ab cd = \[{U + tJ)e~ 2A T (U -U)-{U- U)e 2A T (U + U)\ ab cd , (2.17c) 

M abcd = 1 [([/■ + [7) e -2A T(JJ + - tJ)e 2A T {U - U)\ abcd . (2.17d) 

Then, if we can find the Autonne-Takagi factorization (2.9), the matrix Mmn can be obtained 
in the above manner. Some specific examples are provided in the next sections. 


2.2 Potential 

Having illustrated the way to find out the scalar matrix Mmn, we shall next provide an 
explicit expression of the potential in the SL(8) gauge-field frame. Due to the property 
M mn = (£IM t £1) mn , we have 


m mn = ( (M~ 1 ) abcd {M~ 1 )° b cd \ ( M abcd -M ab cd \ 

V {M-'Uc d ) \-M ab cd M abcd ) 


(2.18) 


Using this relation, some elementary calculations reveal that the potential (2.3) is simplified 
in the SL(8) gauge-field frame to 


2 r 


V = 


9 _ 
672 L 


M abcd (6 be 9 dg W ea ac + 21 e ad 6 bc ) + M abcd (Z be Z da W ac eg + 21 e d i bc ) 


+M ab cd 


ag 


+ 21 8 a c (0 ■ Z) b d ) + M ab 


cd (n cbe 


3 w a9 ec + 21 5 c a (ti • 9) b d 


(2.19) 


where 


W ab cd = 2(M aeb fM ced f + M ae df Mjf ). 


( 2 . 20 ) 


At the origin of the scalar manifold, we have M abcd = M abcd = 5[ a c S b ] d and M ab cd = M ab cd = 
0, recovering the result [41, 43] 


U 0 = ^ [2tr(6> 2 + £ 2 ) - tr(0) 2 - tr(£) 2 ] . (2.21) 

It follows that the potential can be computed straightforwardly with the formula (2.19), 
provided we can implement the diagonalization of 5? (2.9). 

2.3 Kinetic term 

In terms of the mixed coset representative V = the ?]-selfdual part of the T-tensor 

V^jM can be defined by [38] 

Vrtki = [V" 1 ^ - gA™ X M )V] ljkl . (2.22) 
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This quantity defines the scalar kinetic term through (2.4). The Autonne-Takagi factorization 
(2.9) leads to 


fiijkl ( UK/j, U^ijkl j 


(2.23) 


where 

Kf, =d^A — sinh A T UA s ^ sinh A + cosh AU^ A s cosh A 

- coshAt/ t (A% - d^UUftU sinh A + sinh A T U(A a ^ - d^U T U)U cosh A , (2.24) 

with 

A% = ^ m ,(X m + r X M ), A\ = ^4 m ,(X m - T X M ) . (2.25) 

Here, Xjw is a 28 X 28 matrix X Med. = (X a b c d e ft X abcd e f). Using these expressions, one can 
obtain the scalar kinetic term by (2.4) without any difficulty. 


2.4 Gauge kinetic functions 

Next, we proceed to find expressions for the gauge kinetic functions in terms of matrices A 
and U. Let us first focus on the case of the electric gauging, for which the Lagrangian for the 
vector fields is given by [16, 37, 38] 

e~ l C F = -hie(iN A ^)F A ^F^ - \lm(iN A ^)F A ^ * F^ , (2.26) 

where (*T a )^i/ = (1/2 )e flupa F Ap<T and A, X = 1,..., 28 denotes the SL(8) antisymmetric index 
pairs [ab]. In terms of the mixed coset representative 

V m -=(5“ 1 )m-^p-, (2.27) 

the gauge kinetic function AAe = -/V( AS \ is defined by 

V e<j 'JVae = -V A ij . (2.28) 


In the unitary gauge (2.6), we immediately obtain 

r ^ (i + ft tjh (ftftft)^ (i - ft kiMvWo) r cd 

where 


iX abcd = tr 

16 


(2.29) 


tjh(,z) = z 1 tanh(z). (2.30) 

Equation (2.9) yields 

iN = (1 + U tanhAt/ t )(l - UtanhAt/ t ) _1 , (2.31) 

Since the hermitian matrix 4>ft is positive and ^(^^i) -1 / 2 is unitary, we have 

| tanh(\/ 4>(j)^)\ < 1. (2.32) 
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Suppose that det[l —(ft tjh(y^0t)] = 0 is satisfied, then there exists an eigenvector v satisfying 

tjh(\/ <f)(ft)v = v . (2.33) 

This leads to the contradiction since 

v^v = v I tanh 2 (v/00t)n < v^v . (2.34) 

It follows that the gauge kinetic function IVae is always bounded locally. 

Let us next discuss the dyonic gaugings. In this case, we need an extra Stueckelberg 
tensorial held B afll , to compensate the degrees of freedom for the magnetic gauge poten¬ 
tials [38]. This held can be eliminated by going back to the electric frame by the symplectic 
rotation [63]. Let 


L M N = 


_ ( Lt? T A e \ 
- ^ e ) 


(2.35) 


denote the desired Sp(56,M) matrix for which Q'm 0 = Lm N @n° = T (©a“iO). Under the 
symplectic rotation, the gauge kinetic function IVae transforms into 


^ae = (LA^Nzr ~ -LAr)(-L S r - -L Sn Anr) 1 
In the SL(8)-frame (2.1), the no magnetic-charge condition leads to 

L abd be d[e 6 f] V-L ab d[e ^6 f] V=0. 

Assuming #£ = c ^ 0, this equation can be solved to give 

Tab _ I Tabefn r\ 

A cd ~ - VecVdf ■ 

C 

Combined with the symplectic condition L£l T L = U, we have 
L = 


L- l = 


_ f S ■ uj — c 1 { ( t uj) 1 + S ■ uj ■ (9 A 9 )} 
cuj —uj ■ (6 Ad) 

-{9 A6)- t uj ^{(w)- 1 + {0 A0)- t uj- 5} 

t uj • S 


T 

—CUJ 


(2.36) 

(2.37) 

(2.38) 

(2.39a) 

(2.39b) 


where and Sat, = 'S'(AE) are arbitrary 28 X 28 matrices. Here and in the following, we 
will use the notation (S^ A Seabed = *5(1) \a C >5(2)6]^ for the 8x8 tensors <S7 1)2 ). If we make 
the choices, 


uj = t uj = — (9 A 9) 1 , S = — ( t uj) 1 -(9a9) 1 ■ uj 1 = uj 1 
the gauge helds in the electric frame is 


s? = A’ ab X' ab = 


(A 7 -9) A I 8 


— (A 7 • 0) A I 8 


The gauge kinetic function now reduces to 

N' = N[hs + c(9 a9)~ 1 N}~ 1 


(2.40) 


(2.41) 


(2.42) 


We will use this equation later when computing gauge kinetic functions in the dyonic gaugings. 
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3 S0(3) x S0(3) invariant sector of S0(4,4) gaugings 


Although the prescription developed in the previous section can be used in principle to obtain 
the potential for all seventy scalars with any types of consistent gaugings, this is not so 
practical as it stands. In particular, the most difficult and laborious task in this program is 
how to obtain the matrix U in the Autonne-Takagi factorization (2.9), since it admits a large 
extent of non-uniqueness. In this section, we will restrict the consideration to a subsector 
of -E? 7 ( 7 )/SU( 8 ) coset space in such a way that one can circumvent this complexity. We will 
illustrate this by the SO(3) x SO(3)-invariant sector of SO(4,4) gaugings in this section. Still, 
this system is quite rich and yields physically interesting dynamics of scalar fields. 

In the SL( 8 ) gauge-field frame of the SO(4,4) gauging, the embedding tensor takes the 

form 

0 = 14 © HU) , f = epu © (- 14 )], (3.1) 

Here c is the deformation parameter. According to the analysis given in [36, 46], c can be put 
into the range c E [0,1] due to the automorphisms of the theory. 

Among the 70 generators acting translationally on the scalar manifold, the following six 


exhaust the SO(3) x SO(3) 

invariant generators: 


(1,1) 

: 91 — Si a i 8 ; 

Co 

II 

1 

1=1 

4^ 

(3.2a) 

(1,1) 

: 92 — ^(ei234 + 65678) , 

(3.2b) 

(9,1) 

: 93 = S3 A I 8 ; 

S3 = (I3, —3, Q 4 ) , 

(3.2c) 

(1,9) 

: 94 = £4 AI 8 ; 

Co 

II 

© 

1=1 

co 

1 

CO 

(3.2d) 

(4,4) 

: 95 = S 5 A I 8 ; 

^5 — e4 0 e 8 + e 8 0 e4 , 

(3.2e) 

(4,4) 

: 96 = ^(ei238 + 64567) , 

(3.2f) 


where the boldface letters were intended to denote representations of residual gauge symmetry 
SO(4) xSO(4), and we have used the abbreviation e a b c d = e a Ae;,Ae c Aerf. Here, it is understood 
that each g j is used as (frijki in (2.6) through (2.13) to define a generator in the Lie algebra. 
Our generators are related to those of Dall’Agata-Inverso [46] by 

9 f 1 = 92 , 9 f 1 = 92 , 93 l = x (91 + 93 ~ 94) , g2 l = 95 , 9 p 1 = 93 ~ 94, 9™ = 96 ■ 

(3.3) 

In the next two subsections, we shall limit ourselves to further subsectors to illustrate the 
diagonalization method. We derive an analytic expression of the potential for each case and 
discuss the structure of critical points in detail. The full analysis will be done subsequently. 

3.1 SO(4) x SO(4) invariant subsector 

Let us start with the simplest example, the SO(4) X SO(4) invariant subsector, which can be 
treated without using the method explained above and provides an example to be used to 
check the result for a larger sector obtained by the above method. In this sector, only the 
generators g\ and g 2 survive. Let us denote the scalar fields in this subsector as 

(f> ~ xigi + yig 2 , (3.4) 
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where x\ and y\ are the S0(4) X S0(4) singlet scalar and pseudoscalar, respectively. 

Since the matrix cf> can be regarded as a map from M 8 AM 8 ~ M 28 onto itself, it is 
convenient here to decompose the indices as 


M 8 ~ Vi234 © ©5678 , 

R 8 Al 8 ~ [V1234 A 11234] © [^5678 A ©5678] © [^1234 A V5678] ~ M 6 © M 6 © M 16 , (3.5) 


where V\ 234 is the vector space spanned by 1 — 4 indices and similarly for Vs678- It follows 
that the generators g\^ can be decomposed into each block as 

gi ~ Ig © (—He) © ©16 ) 92 — *4 © *4 © ©16 ; (3-6) 

where *4 is the map [* 4 X] ab = \t a bcdX c d onto X ab = X^ aV[ E M 6 , thence 

4> ~ (xi + iyiM)^6 © (—xi + *yi*4)l6 © ©16 , $$ — r\ [16 © 16 © Om]. (3.7) 


Here we have defined x\ + iy\ = ne 1 ® 1 . Instead of computing the Autonne-Takagi factor¬ 
ization (2.10), it is easier in the present case to directly calculate the matrix Mmn by using 
the above matrix representation for </>. First, the coset representative y decomposes as 
5? ~ © 16 © © He © I 32 , where 


= 


( coshri z jhri 
y z* jhri cosh n 


y 2 = ^i(ah —t -xi). 


(3.8) 


with 


Zabcd, — -O $\a ^b] © 


1 . 

— 'lyitabcd • 


(3.9) 


This yields M ~ M\ © Ig © M 2 © Ig © I 32 5 where 


Mi 


/cosh(2ri) + 2xi jh(2n) 2r/ijh(2ri)* 4 \ 

V 2yijh(2ri)* 4 cosh(2ri) - 2xijh(2 n) ) ' 


M 2 = Mi(xi -> — xi). 

(3.10) 


Putting all together, the potential evaluates to 

Q 2 

V = ^-(1 + c 2 )[2 — cosh(2ri)]. (3.11) 

The kinetic term is independent of the deformation parameter and described by the SL(2, M)/SO(2) 
coset space, 


dsr = ——tr (dMdM x ) = dr? + \ sinh 2 (2ri)d$? . (3.12) 

96 4 

It follows that the potential is constant in the angular direction, i.e, 6\ is a modulus, and there 
is no critical point other than the origin. Both of the SO(4) x SO(4) invariant scalars (xi, yi) 
have the mass spectrum —2 (in a unit of cosmological constant) at the origin, and immediately 
fall into the bottomless valley. This behavior is an obvious difficulty in attempting to realize 
inflation around the origin. 
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Let us next compute the gauge coupling function. We first focus upon the case for 
the electric gaugings (c = 0), and subsequently upon the dyonic gaugings by the symplectic 
rotation [63]. By virtue of 4> — (aq + iyiM)^e © (— x\ + iyi* 4 )l 6 © ©i6> we have 


l±4> [ tjh (\fcjxft) = [1 ± (x — \y-k 4 ) tjhri] © I 6 © [1 =F (aq + iyiM) tjhri] © He © Ii6 

cosh 2 n 


[1 + (=f®i + iyiM) tjhri] 1 = 


cosh(2ri) (xi/ri) sinh(2ri) 


[1 + (=Fxr - WiM) tjhri] • 


Hence 


l + tjh(\/ 4>4> t) l — tjh(\/ 4>4> t) 

1 - i(yi/ri)sinh(2n)* 4 


i -i 


/- 


/+ 


where 


Xl 


f± =cosh(2ri) ± — sinh(2ri) 


It follows that 


e 1 C F = - -tr(i ? i 2 • F 12 ) - 


w-tr(i ? \ • Fi) + ^tr(F 2 • F 2 ) 

J- J + 


+ sinh(2ri) 
4ri 


-^-tr(Fi *Fi) + ^-tr (F 2 * F 2 ) 

J- J+ 


(3.13) 


g >i 6 e 1 ~ ifa/ri ] sinh(2ri) * 4 «i 6 ei 16 , (3.14) 


(3.15) 


(3.16) 


where F\ ab = (1/ 2)e a bcdF A cd (-A = 1, 2) and the alternate tensor e a bcd hr M 6 ~ M 4 A R 4 e -> M 6 
should not be confused with the spacetime volume weight, for which the ordinary Hodge dual 
of the 2-form F^ u is denoted by *F fil/ = ^e (U1 / /0(T F po '. 

The term F\ 2 corresponds to the broken gauge symmetry, which will be massive by 
absorbing the Nambu-Goldstone bosons. Hence, we shall focus on F\ and F 2 , which are 
SO(4) gauge field strengths. Splitting these into the chiral components via the decomposition 
SO(4) ~ SU(2) x SU(2) as 

F A{1) ± = F A 12 ±F A U , F a{2) ± = F a 13 ± F a 42 , F a{3) ± = F A U ± F A 23 , (3.17) 


we find that gauge fields with different chirality decouple as 

3 

F A ab • F A ab = 5> a(0 + • F m + + F A(i f ■ F m ~), 

Z— 1 

3 

F A ab • *F A ab = ^(Ta W + • *F Mi) + - F A q~ • *F A(i) -). (3.18) 

Z=1 

Let us move next to the case of dyonic gaugings. Following the prescription discussed 
in section 2.4, we can eliminate the tensorial subsidiary field by going to the electric frame 
by a symplectic rotation. In this electric frame, we have 


N_ N+ 

N' = - -— © I 6 © V-© I 6 © 

1 — 1 — idN+ 


—1 

1 + ic 


© Ii6 j 


(3.19) 


where 


iN± 


1 — i sin 6 1 sinh( 2 ri )*4 

T± 


(3.20) 
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3.2 2,2 truncation 

Let us next discuss the Z 2 invariant subsector studied in [46], which we denote by £ 2 - In 
this case, the generators except for <73 — <74 and g$ are odd under the reflection, hence the 
surviving scalars are described by 


4 > - ®2(ff3 - 54) + 2/256- (3.21) 

It is therefore advantageous to decompose the indices as 
R 8 — V123 © V4 © V567 © V% 1 

R 8 A R 8 — [(V123 A V123) © (V123 A Vs)] © [(V567 A V567) © (V567 A V4)] 

© [(1123 A V567) © (V123 A V4)] © [(V567 A V 8 ) © (V4 A Vs)] . (3.22) 

The scalar held </> is then expressed as 

cj) ~ (x 2 a 3 + iy 2 0i) © I3 © (x 2 cr 3 - iy 2 0i) © I3 © [%2h, ~%2h, ~ 3 x 2 ] • (3.23) 


Since the point at issue is now reduced to the 2x2 eigenvalue problem, the Autonne-Takagi 
factorization ( 2 . 10 ) can be done easily, 


X2 22/2 

22/2 -x 2 


= U 


X 2 - 2/2 


L u. 




(3.24) 


X2 + 2/2. 

Plugging these results into equations presented in section 2.1, the potential is then given by 


n 2 

V = d_ e ex2+2y2 [(1 + c 2 e~ 12x2 )( 1 + e~ 2y2 ) 2 + 3e~ 4 * 2 {2(l - e~ 2y2 ) 2 (l + c 2 e ~ 4x2 ) 

— (c 2 + e~ 4x2 )(l — 6e _2y2 + e _4y2 )}] . (3.25) 

This recovers the potential in [46] after the substitution X 2 —\\n.x, 7/2 —> — |l nr with 
g 2y/2. The target space metric (2.4) is R 2 and reads 

ds^ = 3dx% + dy\ . (3.26) 


As discussed in [46], the potential (3.25) has an off-center critical point for which the mass 
spectrum hinges on the deformation parameter c, and can be small and negative. 

The gauge kinetic function can be evaluated following the argument given in section 2.4. 
In the electric gaugings, we have 


where 


N ~ N + © I 3 © N- © I 3 © [e 2 * 2 I 9 , e~ 2x2 Is, e~ 6x2 J , 


(— ie 2x2 sech(2?/2) ©tanh(2t/2) \ 

V q=tanh(2y2) — ie _2:r2 sech(2?/2) ) 


(3.27) 

(3.28) 


In the dyonic gaugings, the symplectic rotation to the electric frame (2.36) allows us to recast 
the result into 


N' ~ N' + © I 3 © N'_ © I 3 © [/+(2x 2 )l9, /-(—2aj 2 )l6, /+(—6x 2 )] , (3.29) 

where f±(x) = —ie x /{ 1 ± cie x ) and 

N , = _ 1 _ ( c ^ie 2x2 sech(2y 2 ) ^tanh(2y 2 ) \ ^ 3Q x 

± 1 + c 2 - 2icsinh(2x2)sech(2y 2 ) V + tanh(2y 2 ) -cZe -2 x2 sech(2y 2 ) / 

The gauge kinetic function is bounded in either case. 
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3.3 Full analysis 

After having described the two subsectors in detail, we now come to the analysis of the full 
SO(3) x SO(3)-invariant sector. The decomposition of indices is as follows: 

M 8 — V 123 © V 4 0 V 567 © ^8 ) 

M 8 A M 8 —[(V 123 A V 123 ) © (V 123 A V 4 ) © (V 123 A Vs)] 

© [(V 567 A V 567 ) © (V 567 A Vs) © (V 4 A V 567 )] © [(V 4 A Vs) © (V123 A V 567 )] • (3.31) 
Then the generators gi-6 are reduced to 


9 1 = 


92 = 


I 2 


10 1 


© I 3 © 


©I 3 © 


-I 2 


10 1 


© I 3 © O10 , 


© I 3 © O10 


93+94 = 

93-94 = 


-h 


0-3 


0 


© I 3 © 


©I 3 © 


-I 2 


© I 3 © 


-3 


hj ’ 


-cr 3 


-2 


© I3 © 010 , 


95 = \ 1 J © I 3 © 

2^1 


1 I © I 3 © 010 J 
2 CT 1 

— i' 


96= 0 © I 3 © 0 © I 3 © O 10 • 


, — 1 


Let us parametrize the scalars as 


(3.32a) 

(3.32b) 

(3.32c) 

(3.32d) 

(3.32e) 

(3.32f) 


0 =xigi + y\g 2 + x 2 (g :i + . 94 ) + y 2 g& + m(g 3 - 94 ) + 2.w 2 g $, (3.33) 

which can be decomposed into 


where 


4> = $1 © I 3 © $2 © I 3 © (-3x 2 ) © (X 2 I 9 ), 


X\ + x 2 + W\ iyi iy 2 

$1 = ( iyi x\ — x 2 — w\ 1 V 2 

iy 2 w 2 x 2 + 2 w\ 

—xi +x 2 - wi iyi iy 2 

$2 = | iyi —xi - x 2 + wi —w 2 

iy 2 -w 2 ~x 2 - 2w\ 


(3.34) 


(3.35a) 


(3.35b) 


The kinetic-term metric around the origin O is given by 

ds^\o = = 3{dw\ + dx 2 ) + dw 2 + dx\ + dy\ + dy 2 . (3.36) 

All we need to do is to perform the Autonne-Takagi factorization of ^> 1 , 2 . This is yet 
actually redundant, since these 3x3 matrices are related in a simple way as follows 

( x\ + wi iyi iy 2 \ 

iyi xi-wi w 2 \ , (3.37) 

iy 2 w 2 2wi J 
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with I\ 2 = diag(l, —1, —1). It therefore suffices to compute 


&o = UA T U , U £ U(3), A = diag(/ii,// 2 ,M 3 ) • (3.38) 


In Appendix A, we describe how to implement the Autonne-Takagi factorization of the matrix 
$0 of the form (3.37). The eigenvalues /J> t can be obtained by the procedure laid out in 
appendix A and turn out to be the roots of the cubic equation Q{Hi) = 0 given in (A. 5). In 
the present case, we have 


Q(h) —H 3 ~~ ( x i + V\ + V2 + 3tcf + w 2 )h 

+ 2w\ + wi(-2xj - 2y\ + y\ + w 2 ) + xi{w 2 - y{) + 2w 2 y l y 2 • (3.39) 


One sees that the one of the roots is not independent on account of the constraint JA m = 0. 
The eigenvector Vj (A. 8 ) takes the form 


/Hi + (xi + wi)m + 2wi(xi - wi) -w%\ 

Vi= -{Hi + 2w 1 )y 1 +w 2 y2 ■ (3.40a) 

\ -{Hi + Xi - Wi)y 2 + w 2 yi ) 


As demonstrated in detail in appendix A, the unitary matrix U can be parametrized by SU(2) 
Euler angles Hence we can view (x 2 , Hi, @i) under a constraint JA Hi = 0 as new variables 
describing six scalar degrees of freedom. We give in Eq. (B.l) the expressions of the original 
coordinates {xi,x 2 ,yi,y 2 ,wi,w 2 ) in terms of the new variables (x 2 ,/j,i,6i). 

After some tedious computations, we obtain the expression for the potential, 


g~ 2 V = (e 6x2 + c 2 e~ 6x2 ) [ 2 C 321 cosh(2^i - 2y 2 ) + 2C 3 i 2 cosh(2/r 3 - 2/zi) 

+ 2C31-1 cosh(2/j, 2 — 2^3) + C300] 

+ (e~ 2x2 + c 2 e 2x2 ) [2(7111 cosh(2^i) + 2C U o cosh(2 /j, 2 ) + 2Ci 0 i cosh(2/r 3 ) + C 100 ] 

+ 2 ( e 2x2 + c 2 e~ 2x2 ) [C' 1U cosh(2/n) + <7( 10 cosh(2 / u 2 ) + C' wl cosh(2/r 3 )] , (3.41) 
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where 


C300 = 


C31-1 — 

C321 = 
C312 = 
C100 = 
Cm = 
Cno = 
C101 = 
Cin = 

Cno = 
Cioi = 


[3 + p 4 - q A cos(40 2 )] - ^ [( 1 _ P ) 2 cos(20 2 - 20 3 ) - (1 + p) 2 cos(20 2 + 

[(1 — p) 2 cos(20 2 — 20 3 ) — (1 + p) 2 cos(20 2 + 20 3 )] 2 , 

—< 7 2 (psin 20 2 cos 0 3 + cos 20 2 sin 0 3 ) 2 , 

64 

C 3 2i(26 ) 3 -» 20 3 + 7 r), 

9 

32 ’ 



3 

-—sin 2 0 i cos 2 0 3 , 

64 

3 

-wt sin 2 0 i sin 2 0 3 , 

64 

- cos(20 2 ) , 

[ 2 q 2 cos( 20 2 ) - (1 - p ) 2 cos( 20 2 - 20 3 ) - (1 + p ) 2 cos( 20 2 + 20 3 )] , 

Izo 

C(io(203 —> 2 0 3 + 7r) , 


203 )] 2 , 
(3.42a) 

(3.42b) 

(3.42c) 

(3.42d) 

(3.42e) 

(3.42f) 

(3.42g) 

(3.42h) 

(3.42i) 

(3.42j) 

(3.42k) 


with p = cos 0i and q = sin 0i. Note that we have used the condition /a = 0. At the origin 
of the scalar manifold (pi = x 2 = 0 ), one can get the value Vo = (^ 7 2 / 4 )( 1 + c 2 ), as expected. 

In order to recover the potential discussed in the previous subsections, we need to be 
careful since some of C ^ 1 vanishes in those truncations. Consequently we need to keep track 
of higher-order terms in the computation of V and truncate irrelevant scalars at the final step. 
With this careful procedure we are able to reproduce the correct results for the potential value 
and the mass spectrum at the origin [43]. 

Under the constraint YhiPi = 0> the scalar manifold metric (2.4) is also expressed in 
terms of Euler angles as 4 


1 x ' 

ds\ =3 dxl + - d P 2 i + sin h 2 (//2 - /u)Xi + sinh 2 (p 3 - p i)x| + sinh 2 (/ii - p 2 )xi > ( 3 -43) 

i= 1 


where Xi is an SU(2) invariant form 


Xi = d6 3 + cos0id0 2 , 

X 2 = — sin 0 3 d0i + cos 0 3 sin 0i dd 2 , (3.44) 

X 3 = cos 0 3 d 0 i + sin 0 3 sin 0 i dd 2 ■ 


Expanding at the origin, one recovers (3.36) as desired. In the above coordinate system, the 
target space (3.43) has a manifest invariance M x SU(2). 

4 The six-dimensional scalar manifold AIt (3.43) has a product structure A4t — R x A4s, where A4s is 
an Einstein space of negative curvature. Note that Mt itself is a symmetric space where the R.iemann tensor 
is covariantly constant. 
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3.4 Truncation of the gauge sector 

If a generic component of the SO (4,4) gauge field does not vanish, the dynamics of the 
scalar fields fail to be closed inside the SO(3) X SO(3)-invariant subsector exp(J?%), where 
= (gi,--- , go) is the six-dimensional SO(3) x SO(3)-invariant subspace of the tangent 
space of the scalar manifold at the origin. Let us find the restriction on the gauge field 
configuration for the dynamics to be closed within exp(y%). 

Let A denote an element of and A be the corresponding Usp representation, i.e., 

e A G £' 7 ( 7 ) , A = A^Ay) G pus P (c 7 ( 7 )) • (3.45) 

For a generic transformation generated by A, the mixed coset representative V transforms as 

v = $}y i y e A y)y = ,x ( )e A y = (3.46) 


where G SU(8), <5^, y G 74(64) n Usp(56) with 44(n) being the set of hermitian matrices 
of degree n. Let us express y in terms of the exponential mapping as 

y = e x = 1 + X + i X 2 + • • • , X G y (44(64) n Usp(56)). (3.47) 

It follows that the symmetric part of the matrix e A y can be expanded as 

{ e ^y) sym = 1 + X + A + + — [A_, X] + { X , A + } + — A^. + • • • , (3.48) 


where A± = (A ± At)/2. Keeping the terms that are leading order with respect to X and 
linear in A±, the tangent vector X at the origin transforms under the infinitesimal £ 7 ( 7 ) 
transformation A as 


SX = A + + ^[A_,X}. 
In particular, for e A G SL( 8 ), we have 


A = 

A + : 


A A 1 S A 1\ 
S A 1 AM] 

f 0 S A1 

V5A1 0 


G A4(56), 


A = - t A, S 

A A 1 0 \ 

0 A A l) ' 


T S G M(8,R), 


(3.49) 


(3.50) 


For A G so(4,4), the 8 x 8 matrices A and B take the following 4x4 forms 


A = 



S 



(3.51) 


where A\^ are anti-symmetric matrices and B is an arbitrary 4x4 matrix. 

Specializing this transformation formula to the X = 0 case, we find that A+ G is 
required for e A to preserve hence A + oc < 75 . Further, ^6 exhausts all the SO(3) X SO(3)- 
invariant vectors in the tangent space at the origin of the scalar manifold, and SO(3) X SO(3) 
is a maximal proper subgroup of SO(4) x SO(4) containing A_. Hence, the set of all possible 
A_ is identical to so(3)0so(3). Thus, we have found that the gauge fields should be restricted 
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to the sum of the two SO(3) gauge fields and the Abelian gauge held A^ 4,8 ! corresponding the 
generator g§. 

Among these gauge fields, the SO(3) gauge fields do not couple to the scalar fields in the 
sector. Hence, the consistency is trivial for them. We do not consider the Abelian gauge 
held At 4,8 ] corresponding to the Goldstone direction because it has to vanish in the spatially 
homogeneous and isotropic universe. 

Thus, we only have to add the gauge kinetic terms to the Lagrangian for the two SO(3) 
gauge helds 

7l (+)/ = \e IJK A^ JK \ (3.52a) 

A^ 1 = l ei j K A^ +J ^ + ^\ (3.52b) 

where I, J and K run over 1, 2 and 3. The corresponding held strengths can be written 

F (±)/ = dA (±)/ ± | e IJK A {±)J A A (±)A . (3.53) 


The gauge kinetic term is expressed in terms of these helds as 

/e 2 e _1 JSfi = [ Im ( N e)F (t) ■ F (e) ~ Re (l\QR (e) • F (e) 

e=d= 


(3.54) 


From the general formulae in section 2.4, we can find explicit expressions for the gauge 
coupling functions N e as given in Appendix C. 


3.5 Critical points of the potential 

Using the explicit expression of the potential (3.41) in terms of the Euler angles, we now move 
on to finding the critical points of the potential. To this end, the “go to the origin” approach 
by the SL( 8 ) transformation [41, 42] does not work by the following reason. The orbit of SL( 8 ) 
is 35 dimensional, which only covers the half of the coset space F? 7 ( 7 )/SU( 8 ). The overlapping 
orbit of SL( 8 ) and the SO(3) X SO(3) invariant sector consists of the four dimensional space 
corresponding to y\ = y 2 = 0. It follows that the critical point with ( 2 / 1 , 2 / 2 ) 7 ^ (0,0) cannot 
be mapped to the origin by any SL( 8 ) transformations. On account of this fact, we are forced 
to scan the critical points numerically. Unfortunately, the exact critical points we found are 
only the origin and the one found by Dall’Agata-Inverso [46]. 

It is enlightening to derive the mass spectrum at the DI critical point within our formal¬ 
ism, although this has been already addressed in [46]. For this purpose, we define 


y= 2 (^2-^3), z =2 (^2 + ^3). 


It then follows that the DI critical point [46] corresponds to 

7r 2 

01 = -^, 02 = 0 , 0 3 =-ir, x 2 = x *, y = y*, z = 0 . 

where A* = e~ 2x * and Y* = e 2y * satisfy [46] 


/ 1 + 6X2-3 Xf 
X |(3 - 6X| - X 4 ) ’ 


K = 


(-3 hF 2^2)X 4 + 2X2 _ 3 ± 2^2 
X 4 - 6X2 + 1 


(3.55) 


(3.56) 


c = 


(3.57) 






On the Z 2 -invariant plane £ 2 , we introduce the following basis in the tangent space of 
the scalar manifold Tiq = exp(Jt%): 


Id Id Id 

61 V 6 dx 2 ’ 62 \/2 dy' 63 ^6 dz ’ 

15 15 15 

\/2 sinh(2y) 503 5 sinh y 502 ’ 6 y / 2sinhy50i 


(3.58a) 


This basis is orthonormal with respect to the metric K a p of the scalar manifold such that the 
scalar kinetic term reads 


ds 


2 

T 


— K 2 e 1 J£k = ^Kapd(j) a d(j^ , 


(3.59) 


where ((f> a ) = (x 2 , y, z, 63 ,0 2 , 0i). In terms of this basis, we define the 6 x 6 r/ matrix by 


(riab) = (^e a a e b fi D a Dpv\ , (3.60) 

where D a is the covariant derivative in the scalar held manifold with respect to the metric 
K a p. Note that the projection via vielbein is necessary to obtain the mass spectrum for the 
canonically normalized fields. Evaluating at the DI critical point, we find that (r) ab ) has the 
following block-diagonal structure: 


( Jlab )* 


2(3A-j—2A-g+3) 
3(Aj-6A, 2 +l) 



-2 


/ 4(X?-3)(3X?-1) _2_\ 

3(A'^—6A'j + l) ^3 

\ ~73 0 / 


A simple diagonalization gives the eigenvalues 


( MXj-AX^+l) ' 
X*-6X?+1 U 

V 0 0 


(3.61) 


_ -3X 4 + 2Xl ~ 3 ± \J 33A® - 300X6 + 934X 4 - 300XJ + 33 
1,2 “ 3(X 4 - 6*2 + 1) 

2 3X 4 - 10 Xl + 3 ± 2^3X8 - 24X* 6 + 58A , 4 - 24X| + 3 
3,4 “ 3 X 4 - 6X* 2 + 1 ’ 

4(X 4 - 4*2 + 1 ) 

115 x 4 - ex, 2 + 1 ’ 

V 6 = 0. 


(3.62a) 

(3.62b) 

(3.62c) 

(3.62d) 


With the replacement X 2 —> x*, this recovers the result in [46] up to the two different signs 
in front of the square-root. In the critical limit c —> (\J~2 — 1) _ , r\\ and 774 go to zero from 
below, while r \ 2 and 773 positively diverge as shown in Fig. 1. 


4 SO(3) x SO(3) invariant sector of SO(5,3) gaugings 

Following the parallel argument as in the SO(4,4) case, let us next consider the SO(3) x SO(3)- 
invariant sector of the SO(5, 3) gaugings. The embedding tensor is given by [43] 

/prr 

9 = AI 5 © (-l/A)I 3 , t = c9-\ (4.1) 
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Figure 1. The dependence of the ?;-parameters on the deformation parameter c. In the critical limi 
c —> (a/2 — 1)_, r]i and r] 5 are vanishing, whereas others tend to positively diverge. 


The electric gaugings are specified by c = 0, whereas c > 0 for the dyonic case. Note that we 
employed a notation different from that in [43]. The above parametrization (4.1) ensures that 
the origin of the scalar manifold becomes a critical point of the potential, which corresponds 
to the unstable de Sitter critical point. 

In the dyonic case, ref. [36] proposed that the inequivalent gaugings are characterized 
by the eigenvalues of the tensor classifiers 


B at 3 lS = @M 1 ®N S @p e ®Q < *q MNPQ riaerip(; , (4.2) 


where r] a p = tr (t a tp) and qMNPQ = q a ^{ta)(MN(tp)pQ) are the Cartan-Killing metric and 
the quartic invariant of £ 7 ( 7 ), respectively. We follow the convention of the £ 7 ( 7 ) generators 
( ta)M N = {ta)MpQ PN as in ref. [64], Since the exact expression for the eigenvalues of B a / g 7 ' 5 
are not so illuminating, we only comment the results. We found 11 inequivalent eigenvalues, 
all of which are invariant under the inversion c —> 1/c. This renders us to focus on the range 

ce [0,1]. 
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(4.3a) 


The S0(3) X S0(3) invariant generators are given by 

5 8 

(1.1) : 9l = Sr A Ig , Sr = 3 J>“) 2 - 5 ]T(e a ) 2 , 

a= 1 a=6 

3 5 

(14.1) : <7 2 = 5 2 AI 8 , 5 2 = 2^(e“) 2 -3^(e“) 2 , (4.3b) 

a= 1 a=4 

( 14.1) : g 3 = S 3 A I 8 , 5 3 = e 4 0 e 4 - e 5 <g> e 5 , (4.3c) 

(14,1) : g 4 = S 4 A I 8 , S 4 = e 4 <g> e 5 + e 5 <g> e 4 , (4.3d) 

(5.1) : <75 = ^(e 1234 + e 5678 ), (4.3e) 

(5,1): 56 = l(e 1235 -e 4678 ), (4.3f) 


where the left end of each row specihes the representation of SO(5) X SO(3) to which each 
generator belongs in the full sector. In terms of these generators, the scalar fields in the 
invariant sector can be expressed as 

4 2 

= 22 Xigi +22 y^+i ■ ( 4 - 4 ) 

i= 1 j =1 

4.1 Z 2 truncation 

Before going into the full analysis, let us first discuss a simpler case for which some of the 
scalar fields are vanishing. Consider the following index permutations 

(1, 2, 3,4, 5, 6, 7,8) (1, 2, 3,4, -5, -6, -7, -8), (4.5) 

under which the embedding tensor is invariant and the generators 54 and g§ are odd. Then, 
setting X 4 = y 2 = 0 is the consistent truncation. Decomposing the indices as 

M 8 A M 8 ~ [Vf 2 3 A V 423 © Vl23 A V4] © [V(j78 A Vqj 3 © V 5 A V678] 

© [T 123 A V 5 © V 4 A Vg7s © V 4 A V 5 © 4 i 2 3 A Vg 7 8 ], (4.6) 

we have 

(j)~ <g> I 3 © $2 ® I 3 © $3 , (4.7) 

where 


3xi + 2x2 iyi \ 

iyi 3 xi - \x 2 + 5X3 ) ’ 

-5xi iyi \ 


-oxi ig1 

iyi -xi - |x 2 - 4 x 3 


l >3 = [(3xi - 5 X 2 - 2 X 3 )I 3 , (-xi - §x 2 + 2 ^ 3 ) 13 ,3(xi - x 2 ), -(xi - x 2 )H 9 ]. 


(4.8a) 

(4.8b) 

(4.8c) 


In terms of new variable dehned by 

Xi = | (12xi + 3x 2 + x 3 ), X 2 = | (4xi + x 2 - x 3 ), 

X- = x 2 — xi, re ld = X\ + iyi , (4.9) 
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we can see that 


= x-cr 3 + , $2 = x-crz - $q , 

$3 = [(2X 2 - S-)I 3 , (-2X 2 - X_)I 3 , -3x_,x_I g ] , 


with 


$0 


(X x +X 2 vy x \ =JT ( r + X 2 

V iyi Xt-XoJ V -r + X 2 




1 / Vr±X±~ VF=Xj \ 

y/2r \iy/r ~ Xi i\Jr + X x ) ’ 


C/ T C/ = cr 3 . 


The kinetic-term metric is given by 


I 2 x R 2 : 


ds\ = dr 2 + - smh 2 (2r)d6 2 + 3(dX 2 + hx?_), 

15 1 

~ 15dx 2 + — fix! + -dxl + dy\ , 


(4.10) 


(4.11a) 

(4.11b) 


(4.12) 


where at the second line we have Taylor-expanded the metric around the origin. The potential 
is given by 


V = g 2 (V a) + c 2 V {2) ), V {2) = V {1) (c^ l/c,Xi^-Xi,x-^-x-), (4.13) 

where 

I'n, =lTl(e 4 ^ - AV 4 * 2 - A 2 e 2 *"/_ + e~ 2 X 2 U) 

+ ^ (6 - AV 2 **/- - 4 e 22f >/ + ) + + e“ 6Xl / + - 2), (4.14) 

and 


f± = cosh( 2 r) ± cos#sinh( 2 r). 


(4.15) 


The potential for the electric gaugings is therefore given by V = <? 2 V(i)- Expanding the 
potential around the origin, we can get the correct mass spectrum derived in [43] 


V - V 0 ( -30x? + 5xl + - ^yj 


V 0 = 


3 g 2 {c 2 + l)* 


4(3c 2 + l)(c 2 + 3) 


(4.16) 


The existence of critical points which lie outside the SL( 8 )-type critical points is y± 7 ^ 0, i.e, 
6 7 ^ nir (ti £ N). In both of electric and dyonic cases, our numerical scan could not find any 
critical points other than the origin. 


4.2 Full analysis 

We next turn to evaluate the full SO(3) X SO(3) invariant sector. Let us decompose the 
indices as 

R 8 A R 8 ~ [V 123 A V 123 © L 123 A V 4 © V 123 A V 5 ] 

ffi[Vf)78 A Vg 78 © V 5 A Vfi78 © V678 A V 4 ] © [V 4 A V 5 © V 123 A V(j78] , (4.17) 
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leading to 


gi= y 3 J © l 3 

f 2 

92 = -1/2 


1 I 3 © [3, — 19 ] , 


(4.18a) 


I 3 © —3/2 


I 3 © [ 3, 19 ], (4.18b) 


93 = 1/2 


I3 © 1/2 


© I 3 © Qio , 


94 — 0 1/2 I ©I 3 

V 1/20/ 


0 1/2 I © I3 © O10 ■ 

1/2 0 J 


g§= i 0 © I 3 © i 0 © I 3 


t° M 

96 = I 0 I © I 3 


0 i \ 

0 © I 3 © O 10 • 


(4.18c) 


(4.18d) 


(4.18e) 


(4.18f) 


Hence, the coset representative in the invariant sector can be expressed in terms of the six 
scalar fields X{ (* = !,••• ,4) and yj (j = 1, 2) as 


where 


<t> = ^2 Xig.i + ^ Vj94+j - $1 © I 3 © $2 © I 3 © $3 , 
i =1 j =1 


( 3xi + 2x 2 iyi iy 2 

vy\ 3xi - lx 2 + \xz \x 4 
iy2 \x 4 3 xi - lx 2 - 4X3 


/ —0x1 iy 1 11)2 

d> 2 = iyi -xi - |x 2 - 5X3 -lx 4 

\ iy 2 ~2 X '4 -Xl - |x 2 + \ x3 


(4.19) 


$3 = (xi - x 2 ) x [3, -I 9 ] . 
The kinetic term at the origin O reads 


d-s^lo = - [12dx 2 . + 3dx^_ + dx\ + dx\ + 4 (dy\ + dy%) 


where we have defined 


x + = 4xi + X 2 , x_ = X 2 — Xl . 


(4.20) 


(4.21) 


(4.22) 


(4.23) 


(4.24) 


Following the same argument as in the SO(4, 4) case, the 3x3 matrices and d >2 are 
related in a simple manner 


$1 = X_/l,2 + $0 , $2 = X_/l,2 - <J>0 , 


(4.25) 
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with 


/ x+ iyi iV2 \ 

h,2 = diag(l,-1,-1), $ 0 = iyi \(x+ + x 3 ) \xa I . (4.26) 

\ m \x A \{x + -x 3 )J 

Since the matrix 4>o falls into the family (A. 3), we get 

= UA t U , U € U(3), (4.27) 

where the unitary matrix U has exactly the same structure as (A. 9) and the eigenvector takes 
the form 

/ m ? + x + m + \(x\-xl-x\) \ 

Vi= -{fii + \{x+ - x 3 )}yi + \x±y 2 ■ (4.28a) 

\ ~{fM + \ {x + + x 3 )}y 2 + 3 X 42/1 / 

The eigenvalues y t are the roots of the cubic equation Q{yi) = 0, where 

Q{y) =m 3 - ^ i 3x + + x 3 + x l + A (yl + vl)\ M 

+ ^ [ —x + + ^ x 4VlV2 + 2^3( 2/1 — 2 / 2 ) + x +{ a '3 + x 4 — 2(y^ + 2 / 2 )}] ■ (4.29) 

Hence we see yi = 0. Note that the variable x_ decouples from the matrix U, hence also 
from Q{y). By replacing x 2 by x_ in the SO(4,4) kinetic term (3.43), one obtains the full 
kinetic term for SO(5, 3) gaugings, i.e., 

dsT = 3 dx 2 _ + ^^dy 2 + sinh 2 (/r 2 - M3)x? + sinh 2 (/i 3 - y{)xl + sinh 2 (/n - y 2 )xl • (4.30) 

i 

The potential reads 

V = 9 2 (V {1) + c 2 V {2) ), (4.31) 


wiiere 


Vm = 


2 „6a;_ r 


A 2 e 


"W 32 


+ 


+ 


— 1 + cos 2 0i cosh(2 (y 2 — y 3 )) 
sin 2 0i{cos 2 0 3 cosh(2(/r 3 — y±)) + sin 2 0 3 cosh(2(^i — y 2 ))} 

o -2x_ r 1 

--— 6 — —1 |e 2/il cos 2 0i + sin 2 0i(e 2/i2 cos 2 0 3 + e 2fl3 sin 2 0 3 )| 

64 L X z 

A 2 {e~ 2fl1 cos 2 0i + sin 2 0i(e -2At2 cos 2 0 3 + e -2/i3 sin 2 0 3 )} 


64 l 

+ A 2 {e~ 2/il cos 2 0i + sin 2 6i(e~ 2 ^ 2 cos 2 0 3 + e~ 2 ^ 3 sin 2 0 3 )} 

■ 2 n ) + e^ 2/i3 (cos 2 0 3 + cos 2 0i sin 2 0 3 )} 


+ 


3e 2x ~ r 


{e 2mi sin 2 0i + e 2 ^ 2 (sin 2 0 3 + cos 2 0i cos 2 0 3 ] 

— A 2 {e 2/il sin 2 0i + e 2 ^ 2 (sin 2 0 3 + cos 2 0i cos 2 0 3 ) + e 2 ^ 3 (cos 2 0 3 + cos 2 0i sin' 


03)} 


(4.32) 


with 


V( 2 ) = V(i)(x_ ->• —x-,yi ->• —yi, c —>• 1/c). 


(4.33) 
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A careful trace of higher-order terms at the origin (m = x_ = 0), we can recover the desired 
result [43] 


V ~ K 0 


-30xl + 5x% + - (x\ + x\) — - {y'\ + 2/1) 


3g 2 (c 2 + l) 3 
4(3c 2 + l)(c 2 + 3) ' 


(4.34) 


It comes as a surprise to observe that the potential (4.31) is independent of the Euler angle 
02 . This means that the potential possesses the following flat direction 


OV 

002 


0 


0 


- ^47 - 

OX4 Oxs 


0 0 

+ ( 2/1 o-2/2W— 

0 y 2 Oyi 


V = 0. 


(4.35) 


Note that 0/06 2 is a Killing vector of the scalar manifold (4.30). It would be intriguing to 
understand if there is a profound reason for this unexpected flat direction of the potential. 

We do not attempt to write down the exact form of the coupling of the gauge field 
to the kinetic term of the scalars, since their expressions are quite long. But this task is 
straightforward following the arguments of the SO(4,4) case. 

As far as we scanned numerically, there appear no critical points other than origin. 


5 Inflation in the SO(4,4) gaugings 

Let us now undertake the investigation whether a realistic inflationary solution can be realized 
in a maximal gauged supergravity. As is clear from the expression for the scalar field potential 
in terms of the eigenvalues A = (Aj) and the unitary matrix U in section 2, the potential 
at large fields is always dominated by the term exp(^W 2rijXj) with integer rij and bounded 
coefficients. The integer vector n = ( rij ) runs over all possible combination a + b + c where 
a, b and c are vectors that have only one non-vanishing component with value ±1. Then, 
for any one-dimensional sector A j = Cj(f) with a canonically normalized kinetic term —d(j) 2 /2, 
the exponent of the maximally growing exponential factor is given by 2 n • A = 6c m |</>| > 
(lOS/SS) 1 / 2 ^ from 1 = ^c?/6 < (70/6)c^, where c rn is the maximum value of |Cj \. Thus, 
the maximal coefficient in the exponents is too large for a power-law inflation to take place 
[65] unless all the coefficients of the exponential factors with such large exponents vanish in 
some special direction specified by the unitary matrix U in (2.9), and at the same time the 
potential is stationary in all directions orthogonal to that special direction. This implies that 
the chaotic and the power-law inflation do not occur in the maximal gauged supergravity 
except for accidental cases. Only an intermediate or small scale inflation is likely to occur. 
This includes the hill-top-type inflation [66] near a maximum point and the racetrack-type 
inflation [67] near a saddle point. If an additional condition that all coefficients of the growing 
exponential factor vanish in some direction is fulfilled, the Starobinsky-type inflation [2] might 
also be possible along an asymptotically de Sitter flat direction. We will see later that such 
an unexpected situation really happens in our system. 

Another obstruction against inflation in the maximal gauged supergravity is the fact 
that the modulus of the r/ parameter representing the curvature scale of the scalar potential 
is generally larger than the order unity, which spoils inflation. In order for the inflation near 
the critical point to last for a sufficiently long time, the inflaton field must start from a point 
very close to the critical point. Since quantum fluctuations of order of the cosmic expansion 
rate H would push the inflaton away from the critical point, such a classical fine tuning 
becomes unphysical if |ry| > 1. Hence the ^-problem continues to reside in the maximal 
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gauged supergravity. This should not be confused with the conventional r/-problem in the 
N = 1 supergravity [12], since their origins are quite different. 

As we see in more detail below, the theory obtained by the SO(4,4) gauging suffers from 
this 7] problem for a generic value of the deformation parameter. In particular, the saddle 
point of the potential with the highest symmetry at the origin cannot provide an inflation since 
the negative eigenvalues of the mass square at the origin are independent of the deformation 
parameter and of the order unity. However, the SO(3) 2 -invariant saddle points may alleviate 
the rj problem if we fine-tune the deformation parameter to be close to the critical values as 
pointed out by Dall’Agata and Inverso [46]. We have found such a saddle point only for the 
SO(4,4) gauging so far. We therefore restrict our consideration to the SO(4,4) gauging. We 
also study the influence of the gauge flux to see the possibility of the chromo-natural type 
inflation [56, 57]. 

In this section, we adopt the absolute units such that m p i = l/n = I/V&ttG = 1 and 
parametrize the deformation in terms of s instead of c, following ref. [43]. These are simply 
related by 


g-^sg, c=l/s 2 , (5.1) 

hence 6 = £ _1 = s(l 4 ,— 14 ). In this case, the critical value of the deformation parameter 
c = \/2 — 1 corresponds to s = s c = \/\f T l + 1. We can concentrate on the range s > s c due 
to the reflection invariance s —> 1 /s with a parity transformation of the scalar manifold. 


5.1 Description of the cosmological system without gauge flux 

The description of the cosmological system for Einstein’s gravity sourced by scalar fields can be 
specified by the spatially homogeneous configuration. We assume that the metric is spatially 
flat and homogeneous. Namely, we shall consider the Friedmann-Lemaitre-Robertson-Walker 
(FLRW) universe, 

ds 2 = —J,V 2 {t)dt 2 + a 2 (t)dx 2 . (5.2) 

where the lapse function ,.■¥ = <yV(t) will be set to unity at the level of equations of motion. 
Hence, the metric is specified by the cosmic scale factor a(t ) or the e-folding variable a(t) 
defined by a = e“. The cosmic expansion rate H(t ) is expressed as 


H = 


a 



(5.3) 


Here and in the following, a dot denotes differentiation with respect to the proper cosmic 
time t. 

The six scalar Helds are also dependent only on time. In the most part of our analysis, 
we shall use the coordinate system 


System-1 : 0 ( i)(t) = (x 2 (t), , ^t), 6i(t), 6 2 (t), 6 3 (t)). (5.4) 

In this parametrization, the system is governed by the following effective action 



-^e 2 "d 2 - JVV + ^ (3x 2 2 + + lYl l e b fc l - /Zj)Xfc) 

i i,j,k 

(5.5) 
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where V is given by (3.41) and the constraint f^i = 0 is understood. Since the above 
effective action is invariant under the transformation /r 2 —> —^ 2 , 1^3 —* —M3? we can restrict 
the coordinate region to ^2 — ^3 > 0. 

The variation of (5.5) with respect to JV and yields the Friedmann equation and 
the equations of motion for the scalars, respectively. We do not write them explicitly, but can 
be easily deduced. Although the final expressions are messy, equations of motion appear to 
take the simplest form in the coordinate system However, this coordinate system is an 

analogue of the polar coordinates and is singular at the axes where two of fij's are pairwise 
equal, as discussed in Appendix B in detail. This coordinate singularity prevents us to solve 
the equations of motion numerically around these axes. This problem can be circumvented 
by transforming to the original tangent space parametrization 

System-0 : 0 (o) = (aq, yi, x 2 , V 2 , m, w 2 ) ■ (5.6) 

Unfortunately, we fail to have an explicit expression of the Lagrangian in terms of this coor¬ 
dinate system. We find that the coordinate system cj)^ = (u, v, z, 9\, # 2 ) defined in eq. (B.9) 
partially overcomes both of these two shortcomings. See Appendix B for this technical prob¬ 
lem. 

5.2 Slow roll parameters 

As we mentioned at the beginning of this section, the maximal gauged supergravity with 
SO(4,4) gauging suffers from the 77 -problenr for a generic value of the deformation parameter s. 
Before discussing this problem, we first have to address the definitions of slow roll parameters 
for a multi-component system. For the system under consideration, the curvature of the 
potential is represented by the matrix r] = (ijab) of rank six defined by (3.60). At an extremal 
point of the potential, this rj matrix is identical to the mass-square matrix of the scalar fields 
normalized by the value of the potential there. To the contrary, its physical meaning is less 
obvious at a generic point. Nevertheless, we will show that it plays an important role in 
the investigation of inflation. In particular, when the inflaton is a single component scalar 
field, 7 ] = rn'^V"/V together with the e parameter defined by e = m^(V'/V) 2 /2 determines 
the duration of inflation and various observational features of the primordial cosmological 
perturbations produced by quantum fluctuations during inflation. In the case of a slow-roll 
inflation for which e, |?y| <C 1 at the time to when the comoving scale corresponding to the 
present horizon scale comes out of the Hubble horizon during inflation, the observed spectral 
index n s of curvature perturbations around the present horizon scale is determined by these 
parameters as 

n s = 1 + 277 — 6e , (5.7) 

where it is understood that the right-hand side is evaluated at t = to- 

For a multiple-component inflation model it is far from obvious how to quantify the 
“slow roll" as in the single inflaton case. We propose below the alternatives corresponding to 
e and 7j in the single filed case. To this aim, let us first define a sequence of parameters 

- - 1 dnH 

n ~ H-+ 1 dt n ' [ } 

The equations of motion arising from the Lagrangian L = (R — 2V) * 1 — Kapdcj) 0 A *dq t boil 
down to 

4> a + + 3H<j)P + D a V = 0 , (5.9) 
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where T 0 ^ is the Christoffel symbol for the scalar manifold with the metric and D a is the 
associate covariant derivative, i.e., D a V = K a ^dpV . Differentiating the Friedmann equation 

n 2 = l + , (5.10) 


with respect to t and eliminating cj) a with the help of the equations of motion, we obtain 
Raychaudhuri’s equation 

H = -h< a ^ a ^. (5.11) 

Let us now introduce the slow roll approximation for which the equations of motion (5.9) 
can be approximated by 

3H<j) a + D a V ~0. (5.12) 

Equation (5.11) yields 

e H = ei = ^ -^K^daYdpV. (5.13) 

Since the right-hand side of this equation coincides with the e parameter in the single inflaton 
case, we define its counterpart in a multi-component system by 

ey = -^K a td a VdpV . (5.14) 

Next, a differentiation of the slow roll approximation of (5.11) with respect to t, we get 

H 2 H ~ -±-K ap d a VdpV • (5.15) 

18 

Eliminating <j> a by (5.12), we obtain 

e 2 - 2ef ~ -2eir] al3 v a v^, (5.16) 


where v a is the gradient vector of the potential normalized to a unit vector with respect to 
the metric K a p\ 


= -D a V/ (l\ a aD a VD^V\ 


1/2 




(5.17) 


Hence, 


W = rial3V a vP ~ £1 


€2 

2ei 


= V H- 


(5.18) 


This suggests that r/y defined in this equation provides a counterpart in the multi-component 
case to the r) parameter in the single component case. 

Observe that it is more appropriate to work with e# and rjjj rather than ey and r/y when 
we estimate the spectral index in (5.7), because the ^-dependence of n s comes from the time 
dependence of the e parameter. So, in the present paper, we use ey and rjy to get a physical 
understanding of the behavior of inflationary trajectories, while we use en and t]h for the 
terms determined only by the time-dependence of H when we estimate n s and r. Of course, 
they should be equal approximately. In fact, in all numerical calculations of the inflation for 
which more than 60 e-folding was realized, we verified that the difference in these values are 
less than 2% (see Table 2 for some numerical examples). 

Another subtlety of a multi-component system is the non-uniqueness of a slow roll tra¬ 
jectory: there are five-dimensional family of slow roll trajectories in the present case. This 
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non-uniqueness may make the observational predictions highly sensitive to the initial condi¬ 
tion. As we will see in the next subsection, this problem can be partially cured as far as 
slow roll trajectories near the Z 2 -invariant 2D plane E 2 are concerned, since there appear 
two attractor trajectories to which the projection of all sufficiently inflationary trajectories 
on this 2D plane converge. In order to understand this behavior more clearly, let us calculate 
the divergence of the unit tangent vector held v a for the slow roll trajectories (5.17): 

D a v a = -^ 7 = (y*“ - . (5.19) 


This equation implies that D a v a becomes minimum in the direction that minimizes v a vPri a p, 
provided that the 77 tensor is constant. This minimum direction is characterized by the 
eigenvector of the 77 tensor corresponding to the minimum eigenvalue. Hence, if the right- 
hand side of the above equation is negative (this is actually the case in our system), one 
expects that each slow roll trajectory converges to the one determined by the condition 

Vapv 13 = VvKapv 13 . (5.20) 

This condition amounts to the following relation 


D a e v = VV y 2eV D a V. (5.21) 

This equation determines the unique trajectory at least locally. We have confirmed that the 
attractor trajectories satisfy this equation with good accuracy for our numerical solutions (see 
Fig. 2). 


5.3 77 problem 

Since the necessary ingredients and equipments are at hand, let us now discuss the 77 problem 
in our system. As is clear from the arguments in the previous subsection, the 77 parameter 
relevant to our multi-component system is the negative eigenvalue of the 77 matrix with the 
maximal modulus. It coincides with 77 y = v a v^tj a /3 on a slow-roll attractor trajectory. We 
will henceforth denote it by r]y at a generic point of the scalar manifold. 

With respect to the orthonormal basis ( 3 . 58 a), the 77 matrix at a generic point of the 
Z 2 -invariant plane T ,2 has a block structure similar to that at the DI saddle point: 



(5.22) 


where each of the three blocks corresponds to the 2D directions tangent to E 2 , the z — 0% 





Figure 2. An example of an inflationary trajectory on £2 for s = s c + 10~ 3 . The green curves 
represent the contours of ey, and the gray and the steel blue curves represent the contours of V with 
V > 14 and V < T4, respectively. The red cross shows the DI saddle point, and the dotted curves 
represent the V = F* contour. The red curve represents the trajectory obtained numerically. We see 
that it first oscillates but soon settles down to a path on which the gradient vectors of ey and V are 
parallel. 


plane and the 8 \ — 82 plane, respectively. The entries of the matrices are 


Mi = 

M 2 = 
A = 
M 3 = 

m 4 = 


9 


32V 

,2 


fi(X) ( 3 MX) -s 2 -\) h(Y) + 2 9l (X)(3 f 2 (X) - 2) 


g 


32V 

V%g 2 

32V [X 

g 2 


g 3 (x) + 

2 x 


6 


- 3s 2 )X + 6s 2 - 


1 

X 


MY), 


w - w MX) 


32V 

,2 


5 " F 


MX) 3g 2 (X) -s 2 - )MY)~ 16 5 i (X) 


g 


32V 


g 3 (x) + ( -9 — 3s 


X+(6 S 2 -|jl 


MY), 


B = — 


V3g 2 


16V 


93(X) + 


s 4 + 2 


X + 


M 5 = 


M 6 = 


9 


64V 

,2 


_ 6 — 3s 4 

^9 3 (X) H- 2 —X + 


2s 4 + 1 
s 2 X 
6s 4 — 3 
s 2 X 


9 


64V 


6 - 3s 4 v 6s 4 - 3 

9*( x ) + ^2 - x + —mt 


(5.23) 

(5.24) 

(5.25) 

(5.26) 

(5.27) 

(5.28) 


MY) - 6 9l (X)(f 2 (X) + 6) 


fi — 8s 4 

MY)-2{g 3 (X) + =—X + 


(5.29) 


6s 4 -3 

~Mx 


(5.30) 
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Figure 3. The left panel is a 3 D plot of the potential on £2 for s = 2 , and the right panel is its contour 
plot. In the left panel, the blue dot shows the location of the SO( 4 ) x SO( 4 )-invariant extremum, and 
the red dots show the location of the DI saddle points. These critical points are shown by crosses in 
the right panel with some examples of trajectories. 



Figure 4. The time evolution of a(t) for the trajectories shown in Fig. 3 . For trajectories arriving 
at the extremum at the origin, a(t) is drawn up to that point. The time variable t is appropriately 
scaled for each trajectory to make the figure clearer. 


where 


fn(u ) = U n + — , 

u n 


X n s 2 

SnPO = + xK 


( 5 . 31 ) 
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Figure 5. The left panel is a 3D plot of r/v on £ 2 , and the right panel shows contour plots of ey 
(green curves) and rjy (= —0.04, —0.06, —0.08, —0.1 for grey curves, = —0.02 for the light blue curve, 
and = 0, 0.02, • • • for yellow curves). Both are for s = 2. Two examples of trajectories starting near 
the DI saddle point (red cross) on £2 are shown by the black and red curves. 


s - S c 

X * 

y* 

F* 

H, 

y 

OO 

-0.191 

0.573 

0.294 

0.313 

[-0.145,2.309; -1.652, -1.633,0.8069; 0] 

10 -1 

-0.304 

1.01 

0.420 

0.618 

[-0.7908, -0.9002; -0.5643,3.37,1.48; 0] 

10" 2 

-0.389 

1.60 

0.536 

0.661 

[-0.3126,8.52; -0.2075, 6.42, 3.66; 0] 

10" 3 

-0.423 

2.19 

0.592 

0.690 

[-0.09968,26.7; -0.05396,24.7,17.4; 0] 

10" 4 

-0.435 

2.77 

0.611 

0.701 

[-0.03130,85.1; -0.01603,83.1,61.3; 0] 

4- 1CT 5 

-0.437 

3.00 

0.615 

0.703 

[-0.01976,134; -0.01003,132,98.6; 0] 

1.4- 10" 5 

-0.438 

3.26 

0.617 

0.705 

[—0.01168,228; -0.005891,226,168; 0] 

10" 5 

-0.438 

3.34 

0.618 

0.705 

[-9.868 • 10" 3 ,270; -4.971 • 10" 3 ,268,200; 0] 

1.4- 10" 6 

-0.440 

3.84 

0.620 

0.706 

[-3.688 • 10" 3 , 722; -1.849 • 10" 3 , 720, 539; 0] 

10 -6 

-0.440 

3.92 

0.620 

0.706 

[-3.117 • 10" 3 ,850; -1.562 • 10" 3 ,853,639; 0] 

h- 1 ‘ 

O 

1 

O 

-0.440 

6.22 

0.621 

0.707 

[-3.115-lO" 5 ,8.55-10 4 ; 

-1.557 • 10~ 5 ,8.55 • 10 4 ,6.41 • 10 4 ; 0] 

0 

-0.440 

OO 

0.621 

0.707 

[—0, + 00 ; —0, Toe, Too; 0] 


Table 1. Saddle point data 


At the SO(4) X SO(4)-invariant critical point at the origin, the T] matrix is simply 

Hence, the eigenvalues become negative in two directions orthogonal to £2 (one in the second 
block and the other in the third block), and positive in the £2 direction. All the eigenvalues 
are of order unity and independent of the deformation parameter s. Off of the center, negative 
directions change at each point, but the modulus of the eigenvalues are of order unity for a 
generic value of the deformation parameter s. 

Fig. 5 shows the behavior of rjy on £2 for s = 2. As shown in Fig. 3, the potential for 
s = 2 has two DI saddle points in addition to the central extremum on £ 2 . The potential 
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rapidly falls off beyond these DI saddle points. Despite this distinguished behavior, nothing 
special happens for the behavior of r]y around these saddle points. There appears a narrow 
strip on which r]y is around —0.02, but outside this strip \r]y\ is large and hence no inflation 
occurs. 

To be precise, there appear two types of trajectories: those falling toward the waterfall 
valley (the red trajectory in Fig. 3), and those rolling down to the central extremum. In 
the former case, the spacetime immediately collapses into the anti-de Sitter phase, while, in 
the latter case, the expanding FLRW universe with A > 0 is achieved during the roll. If the 
inflaton starts exactly from the Z 2 -invariant plane So and the initial velocity is tangent to this 
plane, the inflaton remains inside this surface due to the Z 2 symmetry of the potential, and 
the universe asymptotes to the de Sitter universe. In reality, quantum fluctuations render such 
a fine-tuning of the initial condition impossible and there inevitably exists a deviation from 
the S 2 plane. This deviation becomes substantial as the inflaton approaches the origin, and 
in the end falls suddenly to the waterfall directions perpendicular to £ 2 . Since the potential 
has no lower bounds in our system, it leads to an immediate collapse of the universe into the 
anti-de Sitter phase. When the potential is modified to have a realistic physical minimum by 
the coupling to an appropriate matter sector, one may be tempted to expect that the universe 
would evolve to a hot big-bang universe. However, this hope is not satisfactorily realized since 
the total amount of e-foldings cannot exceed 0(10) for s = 2 (see Fig. 4), as seen from the 
behavior of ey (see Fig. 5) together with the standard formula of the e-folding number N for 
the single-inflaton slow roll inflation, 



It follows that the inflation is unlikely to occur when s is not close to the critical value s = s c - 
The situation changes drastically when we tune the deformation parameter s to be close 
to the critical value s = s c . As discussed in section 3.5, two negative eigenvalues of the 7 / 
matrix go to zero as s —> s c , and the one corresponding to the £2 direction is the smallest 
when s ~ s c . Hence, we can identify this negative eigenvalue as r]y. In the region where 
r)y < 0, trajectories starting from a point close to the £2 plane converge to an attractor slow 
roll trajectory. Hence, one can expect that the case with small \r]y\ at the DI saddle point 
would drive inflation, provided the inflaton starts from a point near the DI saddle point. Our 
numerical computations confirm that this expectation is indeed true. As is clear from Fig. 1 
and Table 1, we see that the tuning of the deformation parameter s to the level s — s c < 10 -5 
is necessary. 

5.4 Dynamics in the Z 2 -invariant 2D plane £2 

Let us describe results of our numerical studies and discuss the condition under which we 
can realize the inflation consistent with observations in the six-dimensional SO(3) X SO (3)- 
invariant sector by tuning the deformation parameter s. As we have already mentioned, we 
cannot construct a cosmological model describing both the inflationary phase and the subse¬ 
quent big-bang phase in the framework of maximal gauged supergravity, since the potential 
does not have a lower bound and there exists no local minimum at which the reheating oc¬ 
curs. In this paper, we only discuss the features specific to the inflationary phase, such as 
the total e-folding number, the spectral index n s of the scalar curvature perturbation and the 
tensor-scalar ratio r. The issue of reheating is an interesting future work. 
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Figure 6. The left panel is a 3D plot of the potential on E 2 for s = s c + 1.4 • 10~ 5 , and the right 
panel is its contour plot. The DI saddle points and the central dS extremum are shown by red and 
blue dots, respectively, as for s = 2. Some examples of trajectories on the E 2 plane are shown on the 
right panel. 




Figure 7. The left panel shows the time evolution of a(t) for trajectories shown in Fig. 6. The 
trajectory and a(t) for the same solution are drawn with the same color. In the right panel, the blue 
inflationary trajectory is shown with the contours of ey (green curves) and rjy (grey, lightblue and 
yellow curves), ey = 1 for the outermost thick green curve and decreases toward inner curves by the 
factor lCP 1 / 2 . The orange and lightblue curves correspond to rjy = —0.01 and —0.02, respectively, 
and the grey curves correspond to rjy = —0.03, —0.04, • • •. The red dot on the trajectory indicates 
the N = 60 point. 
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Figure 8. The left panel shows examples of trajectories with the contours of V (grey curves), ey 
(light green curves) and r)y (= —0.02 (light blue) and = 0 (yellow)) for s = s c + 10 -10 . ey = 1 
for the outermost thick light green curve and decreases toward inner curves by the factor 10 -1 . The 
trajectories with black, orange, blue and read colors start from ( X2,U2 ) = (x*,4), (x*,3), (x*,2) and 
(x*,l), respectively, while the green trajectory starts from (—0.4924,4) where V = 10F*. The red 
dot on the inflationary trajectories marks the N = 60 point again. The right panel shows the time 
evolution of a(t) for trajectories on the left panel. The time variable t is appropriately scaled for each 
trajectory to make the figure clearer. 



Figure 9. These three graphs show the values of rjn , ejj and n s as functions of the e-folding number 
N between the observational point and the end of inflation for the attractor trajectory in the center 
direction. When s is very close to s c , these functions are independent of s. 
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Let us first discuss the case in which the inflaton starts from a point on the ^-invariant 
plane X 2 with a vanishing initial velocity. In this case, its trajectory stays inside E 2 due to 
Z 2 -invariance. The reflection invariance of the y coordinate allows us to confine to the upper 
half of the S 2 plane. We mainly consider trajectories starting around the DI saddle point in 
this upper half plane. 

Figures 6 and 8 depict the behavior of the potential for s = s c + 1.4 x 10 -5 and s = 
s c + 10 -10 , respectively. One sees that the potential for all s > s c has the same qualitative 
structure as that for s = 2 shown in Fig. 3. As s approaches s c , the y-coordinate of the DI 
saddle point increases, and along with it, the waterfall regions turn into very narrow valleys 
with steep side walls. The contours of V passing the DI saddle point consist of two curves: 
one is the x = x* line and the other approaches the one parallel to the former. The separation 
of these two lines decreases in proportion to s — s c . In particular, the DI saddle point runs 
away to infinity at the limit s = s c , and the waterfall region disappears. There remains only 
a narrow infinitely long asymptotically flat aisle whose width goes to zero at infinity. 

The right panel of Fig. 6 shows some examples of trajectories for s = s c + 1.4 X 10 -5 . 
As this figure shows, there appear two attractor slow-roll trajectories as in the case of s = 2. 
Although all trajectories approach either of them eventually, the total inflation rate of each 
trajectory is sensitive to the moment when the inflaton settles down to a slow roll trajectory. 
In fact, as the left panel of Fig. 7 shows, if the inflaton starts from a point with a large 
potential value, it initially oscillates with a large amplitude before it settles down to an 
attractor trajectory. During this oscillatory phase the cosmic expansion is not inflationary, 
because the kinetic energy dominates over the potential energy. As a consequence, the total 
inflation rate becomes small. This places a constraint upon the initial point of the inflaton. 
Nevertheless, this constraint is not so tight than expected. In fact, inflation with a total e- 
folding number larger than 100 is achieved even if the offset of the initial point of the inflaton 
from the DI saddle point is of order 0.1. 

Let us next look at the values of slow roll parameters at the time t = to when the 
comoving scale corresponding to the present horizon scale comes out from the Hubble horizon 
during inflation. From the requirement that the flatness and horizon problems are resolved, 
the e-folding number N from that time to the end of inflation is constrained to be 50 < N < 
70. This wide range of ambiguity comes from the uncertainties of the inflation energy scale 
and the length of preheating phase. In the present paper, we assume that N = 60 corresponds 
to t = to for definiteness. The results for other values of N can be easily read out from the 
data given in the present paper. 

We give in Table 2 the obtained values of the slow roll parameters 77 # and €fj at N = 60 
for five solutions with sufficiently long inflation starting from different points around the DI 
saddle point. The offset of the initial point from the DI saddle point is the order of 77* (the 
value of 77 at the saddle point) with g = 0.01 or g = 0.1. Note that the solutions to the 
equations of motion ((j>(r), a(r ) and 77(r)/77 *) are independent of g , where r = 77* t is a 
dimensionless time variable. Hence the value of g only affects the initial position through the 
value of 77* oc g. When we compare inflationary solutions with observations, the value of 77* 
is determined by the observed amplitude of the scalar perturbation as we explain later. The 
corresponding value of g is around g = 4 X 10~ 6 (see eq. (5.43) below). It follows that the 
deviation of 77* (77 = 0.1) corresponds to 2 X 10 4 times the actual value of 77. 

A very important feature of our results is that t]h at IV = 60 is around —0.02 independent 
of the initial conditions. This is not related to 77*, the value of 77 at the DI saddle point, because 
it is around —0.01 for s = s c + 1.4 X 10~ 5 (see Table 1), hence 77* is half of 77. We can attribute 
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this feature to the fact that all inflationary trajectories immediately settle down to one of 
the two attractor slow-roll trajectories. On each attractor trajectory, the inflation rate is 
determined by the formula (5.33). As a consequence, the N = 60 point is almost uniquely 
determined by ey. Since r/n is almost equal to rjy for a slow roll trajectory, it turns out that 
the value of r] at N = 60 becomes independent of the initial condition. In fact, the values 
of en and 7 ]h at N = 60 are almost equal for the first two solutions converging to the slow 
roll path down the waterfall, and for the subsequent three solutions approaching the center, 
separately. The value 2en ~ 10~ 4 at N = 60 is also consistent with (5.33) (see the right panel 
in Fig. 7). Though, it is rather surprising that the value of n s at N = 60 has a much smaller 
dispersion than t/h and is almost the same for the two different attractor trajectories. What 
happens here is that the variation of e cancels the variation of r] so that n s = 1 + 2rj — 6e 
become nearly constant. 

These interesting features become clearer when the value of s becomes closer to s c . In 
Fig. 8, we display trajectories starting from four different points on the line X 2 = x* with 
y = y o = 1, 2, 3, 4 for s = s c + 10 -10 . We find that for y$ = 1 and 2, the initial point is below 
the N = 60 point on the attractor trajectory with y ~ 2.69, and as a result, inflation lasts 
only for a short period. In the meanwhile, for yo > 3, sufficiently long inflation is realized. 
As is clear from the explanation above, we will obtain the same values for n s and e for any 
trajectories whose initial value of y is between 3 and y*, provided the initial potential height 
is not so high as illustrated by the green trajectory in Fig. 8. This is because the stability 
of the values of slow roll parameters at the ./V = 60 point arises due to the uniqueness of 
the attractor trajectory. In fact, Table 2 shows that the value of n s at the N = 60 point is 
independent of the initial condition with a good accuracy. 

Even if s — s c becomes smaller further, the behavior of trajectories around and after the 
N = 60 point does not change significantly because the DI saddle point runs away further 
and further from that region, and the potential around the N = 60 point converges to the 
value for s = s c . Fig. 9 shows the ^-dependence of the slow roll parameters t]h and en and 
the spectral index n s in this limiting region of s. 

5.5 Off S 2 trajectories 

Next, we report the results for the trajectories starting from a point outside the Z 2 -invariant 
plane E 2 . As we saw in 3.5, the rj matrix has two eigenvectors with negative eigenvalues at 
the DI saddle point: one is parallel to E 2 and the other is tangent to the {x\,w\) plane (see 
the left panel of Fig. 10). Hence, the inflaton tends to move into the direction off of the E 2 
plane when it starts outside this plane. Nonetheless, the offset from the E 2 plane remains 
small until the inflaton approaches the origin around which all the negative directions become 
perpendicular to E 2 . This is because the eigenvalue in this direction is about one half of that 
in the E 2 direction (see Fig. 1), hence the inflaton acquires a large velocity in the E 2 direction. 

We show illustrative examples of this feature in Fig. 11 and Fig. 12. The solution 
displayed in Fig. 11 admits an initial offset of the inflaton in the positive eigenvalue direction 
of the T] matrix. In this case, the inflaton initially oscillates with a large amplitude passing 
through the E 2 plane orthogonally, but quickly settles down to a trajectory very close to 
the attractor slow-roll trajectory on E 2 . In this stage, there remains a small offset from E 2 , 
which grows rapidly at the final stage. For the solution in Fig. 12, the initial offset of the 
inflaton is of order 10 4 i7*(g = 4 x 10 6 ) with respect to the scalar manifold metric K a p and 
in the negative eigenvalue direction of the rj matrix. In this case, the offset from the E 2 plane 
continues to grow from the beginning. As a consequence, the offset becomes of order unity 
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Figure 10. The left panel is a contour plot of the potential on the (xi,wi) passing through the DI 
saddle point orthogonally to S 2 for s = s c + 10 -10 . The right panel shows a 3D plot of the potential 
on the (£ 2 , 2 / 2 ) plane with X\ = — 2,w± = —0.1. 




Figure 11. The left panel shows the 3D trajectory of the solution for s = s c + 10 -10 with the initial 
condition: x = x*, y = A, z = H*(g = 0.006), Q\ = n/2, 0 2 = 0, 6*3 = 37 t/ 4 . The deep blue dot and the 
lightblue dot on the trajectory show the extremum at the origin and the TV = 60 point, respectively. 
The right panel shows the behavior of the six coordinates around the end of inflation. 


before it comes close to the origin. The inflaton thus feels a potential force in the (x 2 ,j/ 2 ) 
direction quite different from that on £ 2 , as illustrated by the right panel of Fig. 10. It follows 
that the trajectory abruptly changes its direction when the offset grows to some critical value 
as shown in Fig. 12. The inflationary stage is thus destroyed by this property. 

This violent offset behavior, however, occurs only at the final stage of inflation and does 
not affect the evolution of inflaton around and before the TV = 60 point (marked by a lightblue 
dot in Figs. 11 and 12). It follows that the observational parameters take values similar to 
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Figure 12. The left panel shows the 3D trajectory of the solution for s = s c + 10 -10 with the initial 
condition: x = x*, y = 4, z = 0, — y sinh(2i/) cos(20 3 ) = H*(g = O.OGjv^Tr, 9\ = 7 t/2 ,0 2 = 0. The 
meanings of the colored dots are the same as in Fig. 11. The right panel shows the projection of this 
trajectory (red) and that in Fig. 11 on the £2 plane. 


those for trajectories on the £2 plane as shown in Table 2 (the last two rows). 

5.6 Flux contribution 

Having discussed the case without flux, let us next examine the effects of non-vanishing gauge 
flux on the inflationary dynamics. The authors of ref. [56, 57] proposed a new mechanism for 
slowing down the inflaton motion. This class of models is referred to as the chromo-natural 
inflation and caused by the non-Abelian gauge held with a Chern-Simons coupling of the form 

e-^cs = ^Tr(F * F ), (5.34) 

Ja 

where 4 > is an inflaton held, and f a is analogous to the axion decay constant and specihes 
the scale where it comes into play. Fascinating features of this model are that it does not 
demand a hat potential and gives quite different observational predictions [58]. The only 
price to pay is that the coefficient A of the Chern-Simons term must be quite large. To be 
specihc, it must be as big as 200 in the Planck units to achieve the 60 e-foldings. In the 
context of extended supergravities, the Chern-Simons coupling of the type (5.34) is naturally 
incorporated, although the effective coupling constant Re (N) depends on the scalar helds in a 
more intricate fashion. It is therefore interesting to see whether a large effective Chern-Simons 
coupling of this sort appears in our SO(4,4) gauging model, in which there are no controllable 
parameters except for the gauge coupling constant g and the deformation parameter s. 

Another motivation to consider non-vanishing gauge flux is its relevance to anisotropic 
inflation models [59, 60]. Anisotropic inflation is a mechanism driven by the vector or 
higher-rank tensor helds which admit a background configuration with an anisotropic energy- 
momentum tensor. Anisotropic inflation is interesting also from the general relativistic point 
of view, since it deserves a novel counterexample to the cosmic no hair conjecture. In general, 
such an anisotropy quickly decays dynamically with cosmic expansion. For example, the free 
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magnetic field B decays in proportional to 1/a 2 due to the flux conservation, where a is the 
cosmic scale factor. A possible mechanism to prevent this kind of decay of anisotropy is to 
allow an inflaton-dependent gauge coupling of the form 

er 1 &x = -\mF-F. (5.35) 

The Maxwell equation reads V = 0, implying that the combination a? f((p)B is 
conserved. Hence a constant magnetic field is maintained if a? f((p) stays constant. Since 
/ is proportional to the inverse square of the gauge coupling constant, this means that the 
inflaton should evolve toward a strong coupling region with time in order for the gauge field 
to drive inflation. It is thus interesting to see if the N = 8 supergravity displays the behavior 
of this kind. 

Let us now examine these possibilities. The argument discussed in 3.4 demands us to 
truncate the gauge sector to two SO(3) gauge fields A(I = 1,2,3). In order to see the 
possibility of the chromo-natural type inflation, we require that the gauge flux is spatially 
homogenous and isotropic. Then, A^ 1 should be gauge equivalent to 

A ±: = a(t)ip±(t)dx I . (5.36) 


By inserting this into (3.54), we find that the gauge contribution to the action is given by 


Si 


/ 


dta 3 

e=± 


^Im(-Ay 


1 



+ |MOl e(N e )ip 3 


(5.37) 


The explicit expressions for Im(iV±) and Re (N±) are given in Appendix C. By virtue of 
the intricacy of these expressions, it seems less obvious to extract physical information. Our 
numerical computations nevertheless give rise to a simple and surprising result. Figures 13 
and 14 show 3D plots and contour plots of Im (N + ) and Re (N+) on £ 2 , respectively. Fig. 13 
implies that the coupling function Im (N+) is constant with good accuracy except for crossing 
strips parallel to 1/2 = ±X2 • Away from these strips, the function Im (N + ) exponentially damps 
to zero as |rc 2 1 and 1 2 / 2 1 deviate from these strips. This implies that those flattened regions 
correspond to the strong coupling regime. The DI saddle point is inside this region. In 
contrast, Im(A^+) becomes of order unity on the strips around 1/2 = ±X 2 ( X2 > 0). These 
strips contain the critical point at the origin. This implies that for inflationary trajectories 
that terminate at the origin, the inflaton must undergo the transition from a strong coupling 
to a weak coupling. It turns out that anisotropic inflation does not work along that kind 
of trajectory, since the strong coupling is needed so as not to dilute the gauge fields. For 
inflationary trajectories that end in the waterfall valley, the gauge coupling function decreases, 
but its total change is not substantial. 

Let us next discuss the behavior of the coefficient of the Chern-Simons term. Fig. 14 
indicates that there exist several different regions where Re (N + ) remains approximately con¬ 
stant in each domain. We see that the higher value of Re (N + ) is at most of order unity. This 
makes it difficult to drive chromo-natural inflation. We confirmed this statement more rigidly 
by solving the equations of motion of the 6-dimensional scalar system coupled with gravity 
and a single SO(3) gauge held. The inflaton trajectory and a(t) for s = s c + 1.4 x 10~ 5 are 
depicted in Fig. 15. This solution starts from the £2 plane with zero velocity. We take the 
gauge fields to have a vanishing initial electric flux (E oc ip + Hip) but turn on a large ampli¬ 
tude of magnetic flux (ip = ip + = 10, ip- = 0). As we can see from the left panel of Fig. 15, 
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Figure 13. The left panel is a 3D plot of Im (N^) on X 2 for s = s c + 1.4 • 10 5 , and the right panel 
is its contour plot. The potential critical points are marked by crosses. 



Figure 14. The left panel is a 3D plot of Re (iV( + )) on X 2 for s = s c + 1.4 • 10 5 , and the right panel 
is its contour plot. 


the inflaton initially oscillates with a large amplitude due to the magnetic flux, even though 
the initial point is close to the attractor slow-roll trajectory. The oscillation quickly damps 
and the inflaton settles down to the inflationary attractor trajectory since the magnetic flux 
decays in a short time (see Fig. 16). A notable difference from the case without flux is that 
the trajectory exhibits behavior similar to the ones with initial conditions outside X 2 plane 
and the total e-folding number is smaller than the ones on X 2 for a similar starting point, as 
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Figure 15. The left panel shows the £ 2 projection of the trajectory of the solution with the initial 
conditions: x = x*,y = y* — H*(g = 0.01) on £ 2 with flux ip = 10 ,ip/H* = —10. The right panel 
shows a(t) of this solution. 



Figure 16. Behavior of the magnetic flux for the solution in Fig. 15. 


shown in Fig. 15. As Fig. 17 reveals, this behavior arises because the oscillatory motion at 
the beginning is not confined in the £ 2 plane. The magnetic flux effectively produces a large 
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Figure 17. The left panel shows the 3D trajectory of the solution in Fig. 15. The right panel shows 
the behavior of the six coordinate <j), 0 s of this solution. 


offset of the inflationary trajectory from the £2 plane. Thus, the initial strong magnetic flux 
affects the final stage of inflation significantly. As far as this solution, however, the IV = 60 
point (marked by a lightblue dot in Fig. 17) is on the attractor slow-roll trajectory below 
the oscillatory path. As a result, the slow roll parameters at the IV = 60 point take similar 
values to those for solutions without flux. Refer to Table 3 for details. 

5.7 Comparison with observations 

Finally, let us discuss the observational consistency of inflationary solutions we obtained nu¬ 
merically in the six-dimensional sector of the maximal gauged supergravity with SO(4,4) 
gauging. First of all, as is clear from the explanations above and Tables 2 and 3, the spectral 
index n s at the IV = 60 point takes values around 0.96, irrespective of the initial conditions. 
This is within the ler range n s = 0.9639 ± 0.0047 obtained by Planck for a wide range of 
initial conditions, provided that the deformation parameter s is very close to or equal to the 
critical value s c ■ At this critical value, we have only one tunable parameter g which deter¬ 
mines the height of the potential. It follows that this consistency with observations is rather 
miraculous. The situation is similar to the Starobinsky model [2] that admits also a single 
free parameter in the absolute units and predicts a unique value of n s . The similarity has a 
clear physical ground. In our model, there exists only one slow-roll attractor trajectory on 
the £2 plane in the limit s = s c . This trajectory extends out to infinity and the potential V 
on that trajectory approaches a non-vanishing constant V *. The difference V — 14 decreases 
exponentially with the coordinate y^. As a result, the effective one-dimensional potential has 
the same structure as that for the Starobinsky model in the scalar-tensor theory representa¬ 
tion. The only difference is that our model has waterfall fields and inflation ends as in the 
hybrid inflation model. In some aspects, our model provides a variant of multi-dimensional 
version of the Starobinsky model. 

The values of n s and its running are fixed by the slow-roll attractor trajectory on the 
£2 plane for a given N, if the deformation parameter is close to the critical value s c . We can 
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[. 756e-4] 

-.0204 

[-.0206] 

.958 

1.05 

a:* 

y* - 77* 

(5 = 0.1) 

-,352e-2 

381 

-.4340 

2.679 

.709e-4 

[,718e-4] 

-.0202 

[-.0204] 

.961 

1.4-10" 6 

0 

X * 

II 

£ 

0 

727 

-.4341 

2.685 

.832e-4 

[.842e-4] 

-.0182 

[-.0184] 

.962 

10-ro 

0 

X * 

3 

0 

108 

-.4340 

2.683 

.855e-4 

[.866e-4] 

-.0181 

[-.0183] 

.963 

0 

X* 

4 

0 

753 

-.4341 

2.690 

.833e-4 

[.843e-4] 

-.0179 

[-.0181] 

.963 

0 

-.4924 

4 

0 

511 

-.4341 

2.685 

.848e-4 

[.858e-4] 

-.0180 

[-.0182] 

.963 

6.36e-3 

a:* 

4 

-2.12e-3 

740 

-.4341 

2.690 

.830e-4 

[,840e-4] 

-.0178 

[-.0180] 

.963 

1.26e-3 

a;* 

4 

1.26e-3 

752 

-.4348 

2.743 

.673e-4 

[.680e-4] 

-.0160 

[-.0163] 

.967 


Table 2. Examples of numerical results for models without gauge flux. For all models, the initial 
velocity <j) a and the initial values of y\ and wi are zero. The value of g for the initial value of X 2 and 
2/2 indicates that 77* (<x g ) is estimated for that value of g. 


s - s c 

Initial Condition 

N t 

TV = 60 point 


Xi 

X2 

2/2 

Wl 


X2 

2/2 

en[ey] 

1 h [Vv] 

n s 

1.4-10" 5 

0 

x * 

y* + 77* 

(5 = 0.01) 

0 

71 

-.4344 

2.713 

.600e-4 

f.608e-4] 

-.0189 

[-.0192] 

.961 


Table 3. An exampel of numerical results for the model with gauge flux. 


thus read out the value of the running —dn s /dN from the right graph of Fig. 9 as 

f -9.3 x 10- 4 (TV = 50), 

^ ~ < -6.7 x 1CT 4 (TV = 60), (5.38) 

dN { -4.6 x 1CT 4 (TV = 70). 

These values are consistent with the observational constraint —0.0033 ± 0.0074 obtained by 
Planck temperature data and the CMB lensing [68]. Thus, the scale-dependence of the 
spectral index is very small in our model. 

The spectral index is independent of the height of a potential. In contrast, the amplitudes 
of the curvature perturbations ^^(k) (representing the dimensionless power spectra) and 
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tensor perturbations £Ph(k) (representing the gravitational waves) depend on the potential 
height through the cosmic expansion rate H as 


& c (k) 


1 (H(t k )\ 2 

2e(4) \2irm pi J 


&h(k) 


H(t k ) y 

2vrm p i/ 


(5.39) 


where we have restored the Planck mass and tk denotes the cosmic time when the comoving 
scale corresponding to the conroving wavenumber k comes out of the Hubble horizon. From 
these equations, the famous consistency relation for the tensor-scalar ratio r follows: 


r = 


~&b 


= 16e. 


(5.40) 


We have to use this formula with care because this formula holds only under a set of as¬ 
sumptions that the slow-roll inflation by a single inflaton with the standard kinetic term and 
negligible curvaton contribution. 

In our numerical solutions, the values of e at the N = 60 point has a larger dispersion 
than that of n s , but still remains of the same order. In particular, for inflationary trajectories 
ending around the origin, our model predicts 


0.96 x 10" 3 < r < 1.37 x 10“ 3 . 


(5.41) 


This is about one third of the value in the Starobinsky model and around the border of the 
sensitivity of the next generation B-mode satellite experiment LiteBIRD [69]. 

Assuming that the major part of scalar curvature perturbations is produced by inflaton, 
the CMB normalization by Planck and WMAP [70], 

^ c (/co) ~ 2 x 10“ 9 , (jfeo = 0.002Mpc _1 ) (5.42) 

fixes the value of the gauge coupling constant g as 

g ~ 4 x 1(T 6 . (5.43) 


6 Concluding remarks 

In the present paper, we have addressed the problem whether the N = 8 gauged supergravity 
with a one-parameter deformation realizes a realistic inflationary cosmology. Our work is mo¬ 
tivated in part by the theoretical restriction and the ultraviolet control of IV = 8 supergravity. 
If cosmological predictions of such a theory that is constrained by general physical principles 
are consistent with observations, it will provide us with a hint for deep understanding of the 
fundamental laws of Nature in contrast to the phenomenological framework provided by the 
N = 1 supergravity. From this perspective, we have made elaborated studies of the dynami¬ 
cal structure in the SO(3) x SO(3) invariant sector of the SO(4,4) and SO(5,3) gaugings of 
N = 8 supergravity, extending the work of [46]. 

In the SL(8) gauge-field frame, we have outlined a general recipe which enables us to 
compute the potential analytically. This strategy is useful for searching “off-center” critical 
points. Although the final expression is quite messy for the full SO(3) X SO(3)-invariant 
sector where six nonvanishing scalars survive, our work may serve as a cornerstone to the 
analytic scan of the critical points on the whole scalar manifolds in extended supergravities. 
Our formulation indeed works also for the potential of the SO(8) gaugings obtained in the 
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literature [36, 48, 49]. The present analytic way of obtaining the potential is applicable as far 
as the scalar fields take values in the coset space, e.g, N > 3 supergravities. The application 
of the present study for these cases is a plausible future work. 

Unfortunately, our numerical computations show that there appear no critical points 
except for the known ones for the SO(3) x SO(3) invariant sector of SO(4,4) or SO(5,3) 
gaugings. In order to find out the other critical points, we need to consider a larger scalar 
sector that does not respect the SO(3) X SO(3) invariance, or other gauge groups. At present, 
no de Sitter critical points have been found in a gauge group other than SO(4,4) and SO(5,3). 
As far as the authors know, however, there seems no no-go theorem which prohibits the de 
Sitter solution in other gaugings. 

Although metastable de Sitter critical points are of importance in the context of dark 
energy, the nonexistence of these extrema is not a serious problem in the context of infla¬ 
tionary model-building, because otherwise the inflation never ends. In light of this, we have 
extensively investigated numerically the dynamics of six scalar fields coupled with gravity 
and gauge fields in the SO(4,4) gaugings. We concluded that a sufficient duration of inflation 
can be easily achieved, provided the deformation parameter is very close to or equal to the 
critical value. The initial point of the inflaton should be around the saddle point found by 
Dall’Agata-Inverso [46], but the offset can be order unity in the Planck units. We have fur¬ 
ther determined the values of the slow roll parameters at the N = 60 point. We have found 
that the spectral index n s for the scalar curvature perturbation is miraculously insensitive 
to the initial condition and the value of initial flux, and in the ler range reported by Planck 
satellite [1]. We have explained why such a miraculous situation arises on the basis of the 
existence of special attractor slow roll trajectories. We have pointed out that the situation is 
analogous to the Starobinsky model and that our model is a multi-dimensional extension of 
the Starobinsky model. In fact, our model predicts a small tensor-scalar ratio around 10 -3 , 
which is close to that of the Starobinsky model. 

We have also explored the possibility of the chromo-natural inflation due to the Chern- 
Simons terms and the anisotropic inflation due to the non-trivial inflaton dependence on the 
gauge coupling function. We have denied these possibilities based on the inspection of the 
behavior of the gauge coupling functions and the numerical studies of the evolution of the 
system with non-vanishing flux. Because it is a common structure of gauged supergravities 
with dyonic gauging that the real part of the gauge coupling function appearing in the Chern- 
Simons term is bounded in the electric frame, we can conclude that the chromo-natural 
inflation fails to occur in the framework of IV = 8 gauged supergravity. 
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A Autonne-Takagi factorization 

In this appendix, we briefly outline how to implement the Autonne-Takagi factorization for 
a 3 x 3 complex symmetric matrix: 

<S> = UA T U, A = diag(/ri,/X2,M3), U 6 U(3), (A.l) 
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where yi are real eigenvalues. Eq. (A.l) is tantamount to the following eigenvalue problem 


$«(i) = , Uij = . (A.2) 

The 3x3 complex symmetric matrix <h that we encountered in the body of text takes the 
universal form 


/ a iyi iy 2 \ 

$ = I iyi b w I , (A.3) 

\ * 2/2 w c ) 

where a, 6, c, w, y \, y 2 are real quantities. Decomposing the eigenvector into real and imaginary 
parts as rq = Xj + iY t . we have 

/ M ~a iyi iy 2 \ / X 1 \ / m + a iyi iy 2 \ / 1 Y 1 \ 

I iyi y-i + b w ) I iY 2 ) + I iyi b - m -w 1 I X 2 I = 0. (A.4) 

\ iy 2 w yi + cj \iY 3 ) \ iy 2 -w c- y t ) \X 3 ) 

Let us define Q(y) as the determinant of the first matrix in (A.4), i.e., 

Q(y) =y 3 + {b + c- a)y 2 + {be - a(b + c) - w 2 -y\- yl}y 

+ a(w 2 - be) - cy\ - by 2 + 2wyiy 2 . (A.5) 

In this case, one can easily verify that the determinant of the second matrix in (A.4) yields 
-Q(-y). Together with 4> ( ld = /r 2 l 3 , we have the property 

0 = details - = -Q(y)Q{-y ). (A.6) 

Setting Q{y) = 0 loses no generality and this choice requires Yi = X 2 = X 3 = 0. It 
therefore turns out that the complex eigenvalue problem (A.l) now reduces to the following 
real eigenvalue problem: 

/ a —yi y 2 \ f Xl \ 

Av = yv, A = I —yi -b —w I , v = I Y 2 J . (A.7) 

\ - 2/2 -W -c ) \y 3 ) 

A straightforward computation shows that 

/ x \(i) \ ( y 2 i + {b + c)yi + be - w 2 \ 

Vi= I y 2 (i) I = I ~yi(yi + c) + wy 2 I GW 3 . (A. 8) 

V *3(0 J \ - 2 / 2 {yi + b) + wyi J 

This gives 

/ Clx m c 2 x 1(2) c 3 x m \ 

U = iCiY m iC 2 Y 2{2) iC 3 Y 2{3) , C, = {X 2 {{) + Y 2 {i) + Y 2 ^) 1 / 2 . (A.9) 

\iCiY m iC 2 Y 3 ^ 2 ) iC 3 Y 3 ( 3 } J 

This form of matrix U fulfills 

U T U = diag(l, —1, —1), Ui T U 2 = 0 , U = U\+ iU 2 , (A.10) 
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leading to 


0 = U! + U 2 , O t O = I 3 • 


(A.ll) 


Note that O is an orthogonal matrix appearing in the diagonalization of 3 X 3 real eigenvalue 
problem (A. 7). A suitable choice of the phase allows us to set O E SO(3), so that one can 
employ a parametrization in terms of Euler angles (0 1 , 62 , # 3 ), 


O 


( cos 9\ — sin 9\ cos 9 3 

sin 6 \ cos #2 cos 9 1 cos 6*2 cos 9 3 — sin #2 sin $3 
sin 9\ sin $2 cos 9\ sin $2 cos $3 + cos 02 sin 63 


sin 9\ sin $3 \ 

— cos 9 1 cos #2 sin 6*3 — sin #2 cos 9 3 

— cos 9\ sin 02 sin 0 3 + cos 02 cos 03 / 

(A-12) 


B Coordinate Systems of the six-dimensional subspace 


In this appendix we discuss the coordinate transformations for the target space in the SO(4,4) 
gaugings. The original coordinate system = ( x h x 2 , 2/1, 2 / 2 , u>i, W 2 ) defined in (3.33) is 
regular on the whole space Eg = exp(J^). From eqs. (3.38), (A. 9), (A.ll) and (A.12), 
one finds that these coordinates can be written in terms of the polar coordinates (f)M\ = 
(x 2 ,n 2 ,^ 3 , 01 , 02 , 03 ) as 


Xl 


V 1 
2/2 

W\ 


W 2 


— 2 C l(l C 2)(-^ + c 3) — 9 + c 3 + 2 C 2^ + C 1 S 2C2S3 Cs| /i2 

+ {“^(1 + ^X 1 + s i) - \ + S 3 + ^(1 + c|) - C1S2C2S3C3 j M 3 , 

Si [{-CiC 2 (l + C3) + S2S3C3} /.i 2 + {-CiC 2 (l + S3) - S2S3C3} fl;i] , 

Si [- {ciS 2 (l + C3) + C2S3C3} (J.2 - {ciS 2 (l + S3) - C2S3C3} /1, 3 ] , 

~2 C 1' S 2( c 1 + 1) + 2 — C 2 + 9 C 2 C 1 — S 2 S 3 ClC 2 C 3 | (12 

+ |-^C?sl(s§ + 1) + ^ - c| + ^c|s| + S2S 3 ClC2C 3 | /13, 

{s 2 c 2 (-Ci(c| + 1) + S3 + 1 ) - C1S 3 C 3 (2C2 - 1 )} 

+ {s2C2( — cf (S3 + 1) + C3 + 1) + CiC3S3(2c| — 1)} /i'j, 


(B.la) 

(B.lb) 

(B.lc) 


(B.ld) 


(B.le) 


where 

d = cos(0j), Sj = sin(0j), (* = 1,2,3). (B.2) 

The right-hand side of this coordinate relation is invariant under the following six transfor¬ 
mations: 


T\ 

: 03 —t 03 + 7T, 


(B.3a) 

t 2 

: 01- 

■> 01 + 7T, 

03 —> -03, 

(B.3b) 

t 3 

: 0x - 

->-0i, 

02 -> 02 + 7T, 

(B.3c) 
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(B.3d) 


7r 

T 4 ■ (M2,M3) (M3,M2), ^3 -> #3 + 7^ 

T5 : H2 —1 — (m2 + M 3 )r 

(si,ci) -s- (±\/l - «?c|, -sic 3 ), 

(«2, C 2 ) Ah(ciS 2 C 3 + S 2 S 3 , C1C2C3 - S2S3), 

(« 3 , c 3 ) ->• -/V 2 (siS 3 , ci), (B.3e) 

2); • M 3 1 —(M 2 + M 3 ), 

(si,ci) -s> (±\/l - s?si,sis 3 ), 

(«2, C 2 ) -> N 1 (-C 1 S 2 S 3 + C 2 C 3 , -(C 1 C 2 S 3 + S 2 C 3 )), 

(•S 3 , C 3 ) N 2 {-Ci, S 1 C 3 ) (B.3f) 

Here, note that the subspaces M 2 = M 3 , Ml = M 2 and M 3 = Ml with Ml + M 2 + M 3 = 0 
are invariant under the transformations T 4 , T 5 and Tq, respectively. Further, each of these 
transformations exchanges two of these three subspaces. Hence, six wedges cut out from the 
M 2 — M 3 plane by these three lines are all equivalent if we neglect the angular coordinates 0j 
(i = 1,2,3). If we select one of the wedges as the fundamental region in the M 2 — M 3 plane 
with the help of the transformations T\, T 2 and T 3 , we can restrict the ranges of the Euler 
angles as 

0 < 01 < 7T, 0 < 02 < 7T, 0 < 03 < 7T. (B.4) 

Instead, if we select two of the wedge regions as the fundamental region, say, the region cut 
out by two lines 2 m 2 + M 3 = 0 and M 2 + 2 m 3 = 0, 

(M 2 ,M 3 )=p( 2 ,-l) + ?(l,-2) = (2p + g,-p-2g), pq > 0, (B.5) 

we can reduce the range of d 3 to 7 r /2 < 0 3 < it with the help of T 4 . In this choice of the 
fundamental region, the SO(3) x SO(3)-invariant plane S 2 can be represented as M 2 + M 3 = 
0, 0i = 7 t/ 2 , d 2 = 0, 03 = 37 t /4 as used in the present paper. 

Because these three subspaces m« = Mj (* > j = 1,2,3) in the coordinate system 
correspond to a single subspace in % and play the role of an axis for the Euler angles, this 
coordinate system is singular at this axis. This can be easily confirmed by calculating the 
Jacobian of the transformation: 

= \ sin( 0 i)(M 2 - M 3 )( 2 M 2 + M 3 XM 2 + 2m3)- (B. 6 ) 

In particular, the x-axis on T, 2 is contained in the singular subspace p 2 = M 3 - 

Although the coordinate singularity at the intersection of these three subspaces are 
quite nasty and can be removed only by going back to the original coordinate system <^ 0 ), 
singularity outside this intersection can be removed by introducing Cartesian-like coordinates. 
As an example, let us consider the fundamental region cut out by two lines M 2 = M 3 aR d 
M 2 + 2 m 3 = 0, and introduce z and r by 

r 

M 2 = z + r, p 3 = z-r, r> 0, 2 > - • (B.7) 

In terms of the variable r, we can introduce the two-dimensional Cartesian coordinates (tt, v) 
by 

u = rcos(29 3 ), v = rsin(29 3 ). (B.8) 
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The original coordinate system <^ 0 ) can thus be expressed in terms of the new coordinate 
system </>( 2 ) = (w, v, z, 9\, # 2 ) as 


w\ + x\ = s\u + (1 — 3 c 2 )z, 

w\ — x\ = {— 2 C 1 C 2 S 2 V + (c 2 c 2 — s 2 )it} + (1 — 3sfc%)z, 
W 2 = ci(l — 2 c|)n — (1 + c\)c 2 S 2 U + 3sfc2S2Z, 

yi = si (s2v - ci02(11 + 3 z )), 

U2 = Si (-C2V - ClS 2 {u + 3 z)). 

The Jacobian of this transformation is given by 

= sin(#i)(9z 2 — r 2 ). 


D <t>( 2) 


(B.9a) 

(B.9b) 

(B.9c) 

(B.9d) 

(B.9e) 


(B.10) 


Hence, the coordinate system 0^ 2 ) is singular only at r = 3 z, i.e., /i 2 + 2/13 = 0 apart from 
the singularity of the polar coordinate at 6\ = 0. 

When we use this coordinate system in the study of trajectories around X 2 . we have 
to map the fundamental region we adopted in the present paper to the above one by the 
transformation T5 or T@. We should be also careful about the fact that this coordinate 
system is still singular on the x-axis on E 2 . When we consider a trajectory that pass through 
this axis, we have to adopt the original coordinate (^(o) near the axis where an approximate 
expression for the action can be obtained by the Taylor expansion. 

The kinetic term of the effective Lagrangian for the scalar fields can be written in terms 
of the coordinate system <^ 2 ) as 


ds\ = dx 2 + du 2 + dv 2 + 3 dz 2 


+ 


jh(2?’) 2 — 1 


(■udv — vduY 


+A]\\{2r) 2 (udv — vdu + r 2 cos(0i)d# 2 ) cos(0i)d0 2 


(cosh(6z) cosh(2r) — 1) (d0\ + sin 2 9\d0 2 ) 

+jh(2r) sinh(6z) {u(d0f — sin 2 9\d9 2 ) — 2xsin0id#ia!# 2 } ■ 


(B.ll) 


The potential for the scalar fields reads 

„2 


g~ 2 V = 


X 3 


+ 


X 3 


D : >i cosh(6z) cosh(2r) + sinh(6z)jh(2r)(T>3 2 'n + D 33 v) 


+jh(2r) 2 (T> 34 it 2 + D 35 uv + D 36 v 2 ) + D 37 


+ 


32 


s 2 X + 


1 


s 2 X 


— c 2 cosh(4z) — s 2 cosh(2z) cosh(2r) 


+ 


— s{ jh(2r) sinh(2z)u + 6 
2X r 3 


X s 
+ 




16 


(— cosh(4z) + cosh(2z) cosh(2r)) s 2 cos(20 2 ) 


+ -jh(2r) sinh(2z) { —(1 + c 2 ) cos(20 2 )u + 2ci sin(20 2 )n} 


(B.12) 
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where 


D 3 i = ^ { s i cos( 46» 2 ) - 1 - c?} , 

D 32 = {l + (1 + c?)sin 2 (2(9 2 )} , 

D 33 = S ^sin(46» 2 ), 

lb 

D 34 = —{cos(4$i) + 28cos(26»i) + 35)cos 2 (20 2 ) 

25b 

— cos(4$i) + 4cos(20i) — 3}, 

D 35 = -^(7cos(6>i)+ cos(36>i))sin(46» 2 ), 

d 36 = 7^7 {(cos(40i) + 28cos(20i) + 35) sin 2 (2(9 2 )} , 

D 37 = ^ (- cos(40i) + 4 cos(20i) - 3) cos(40 2 ) 

+—(cos(40i) + 4cos(20i) + 27). 


(B.13a) 

(B.13b) 

(B.13c) 

(B.13d) 

(B.13e) 

(B.13f) 

(B.13g) 
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C Gauge Coupling Functions for the SO(3) x SO(3) gauge fields in the 
SO(4,4) gauging 

The gauge coupling functions N± in (3.54) are related to the general gauge coupling matrix 
,//' M'l V’A i n the electric frame as 


N + = jY'^ 12 \ 

In terms of 

X = e~ 2x2 , Y = 

N+ can be written as 

Im (N+) = 

with 


X_ = .,y'[56][56]_ 

(C.l) 

e 2 ^ 2 , Z = e 2f13 , 

(C.2) 

Re (N+) = 

(C.3) 


A 


B 


C 


s 4 X 6 Y 2 Z 2 + X 4 (a 442 Y 4 Z 2 + a 433 Y 3 Z 3 + a 424 Y 2 Z 4 + a 421 Y 2 Z 
+a 412 YZ 2 + a 400 ) + X 2 (a 244 Y 4 Z 4 + a 232 Y 3 Z 2 + a 223 Y 2 Z 3 
+ci 220 Y 2 + a 2U YZ + ci 202 Z 2 ) + s 8 Y 2 Z 2 , 

- s 4 x\x 4 {b 423 Y 2 Z 3 + b 432 Y 3 Z 2 + b 4U YZ) 


+X 2 (b 243 Y 4 Z 3 + b 234 Y 3 Z 4 + b 231 Y 3 Z + b 222 Y 2 Z 2 
+b 2 i 3 YZ 3 + b 2 ioY + 6201 -Z) 

+b 033 Y 3 Z 3 + b 012 YZ 2 + b 021 Y 2 z], 


X 4 (c 442 Y 4 Z 2 + c 433 Y 3 Z 3 + c 424 Y 2 Z 4 


+C 4 2 lO 2 Z + C 4 i 2 T Z 2 + c 4 oo) 
+X 2 (c 244 Y 4 Z 4 + c 232 Y 3 Z 2 + c 223 Y 2 Z 3 

+c 220 Y 2 + c 211 YZ + c 202 Z 2 ) + s s Y 2 Z 2 


(C.4a) 


(C.4b) 


(C.4c) 


Here, a^, bijk and are functions of the Euler angle 9i (i = 1,2,3). iV_ is easily obtained 
from N + by 

IV_ = N+(Y -> l/Y, Z -> 1 /Z). (C.5) 

The coefficient functions a^, bijk and are expressed in terms of 

Ci = cos (9i), Si = sm(9i) ( i = 1, 2,3), (C.6a) 


as follows: 


a 442 ={(— 4cfs2 + 4c\s\ — 24c 2 S2 — c\ + 24 c\s\ — 4s\ — 2 cf + 4s\ — 1)03 
+ 8ciC 2S2S3(2C2 — l)(c? + l)c| 

+ (24c?s^ - 24 c\sl + + 2c? - 8 s\ + 2 )c§ 

- 8 cic 2 s 2 S 3 ( 2 c? - l)c 3 - ( 2 s\ - l) 2 }s 8 

+ { — 2(2 S? — l)s?(cj + l)Cg — %Cis\s 2 C 2 S 3 4 + 2c|(2s? — l)s?}s 4 

- 44 (C.7a) 
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0433 = {(8cf4 ~~ 8 cfs| + 48 cfsf + 2 cf — 48c?s? + 8 s? + 4 cf — 8 s? + 2)03 
-16c 1 c 2 s 2 s 3 (2c$ - 1 )(cf + l)c| 

+(— 8 c?s? + 8 c 4 s? — 48c?s? — 2 cf + 48c?s? — 8 s? — 4c? + 8 s? — 2 )c| 

+8cic 2 s 2 s 3 (2c? - l)(c? + l)c 3 - 8 c?s?c?}s 8 
+{4(2s? — l)s?(c? + l) c | + 16ci-sJ-s 2 c 2 s 3 c| 

+(-4s?(2s? - l))(c? + l)c| - 8 cis?s 2 C 2 S 3 c 3 }s 4 

+2c 3 6'i — 2s?c§ (C.7b) 

0424 = {(—4c?s? + 4c?s? — 24c?s? — cf + 24c?s? — 4s? — 2c? + 4s? — l)c| 

+8cic 2 s 2 s 3 (2c? - 1 )(c? + l)c| 

+2c?(4c?s? — 4c?s? + 12s? + cf — 12s? + l)c§ 

-8s 2 c 2 c?s 3 (2c? - l)c 3 - c?(2s? - l) 2 }s 8 
+{—2(2s| - 1 )s?(c? + 1)4 - 8s?cis 2 c 2 s 3 c 3 
+2(2s? - l)s?(2c? + 1)4 + 8ciS?S 2 C 2 S 3 C 3 - 2c?s?(2s? - l)}s 4 
-sfsj, (C.7c) 

0421 = s?{(-8c?s? + 8c?s? - 8s? - 2c i + 8s?)c| + 8 cic 2 s 2 s 3 (2c? - l)cg - 8s|c|}s 8 

+s?{4c?(2s? - l)c| + 8 s 2 c 2 c 3 cis 3 }s 4 - 2s?c?c|, (C.7d) 

0412 = S?{( 8 C?S? ~ 8 c?S? + 8 S? + 2c l — 8 s?)c? — 8 ciC 2 S 2 S 3 ( 2 c? — 1 )C 3 — 2c?(2s? — l) 2 }s 8 
+s?{-4c?(2s? - 1)4 - 8 S 2 C 2 C 3 C 1 S 3 + 4c?(2s? - l)}s 4 - 2 c?s?s§, (C.7e) 

0400 = — (2s? — l ) 25 ! 58 — 2s?c?(2s? — l)s 4 — cf, (C.7f) 

a 2 44 = —cfs 12 — 2 s?c?( 2 s? — l)s 8 — (2s 2 — l) 2 s?s 4 , (C.7g) 

0232 = “2s?C?S?S 12 + S 2 { — 8 s 2 C 2 C 3 ClS 3 + 4c?s?(2s? - l)}s 8 
+S?{( — 8 c?S? + 8 C?S? — 8 s? — 2c l + 8 s?)s§ 

-8cic 2 s 2 s 3 (2c? - 1 )c 3 - 8s?c?}s 4 , (C.7h) 

0223 = -2s?c?c|s 12 + s?{-4c?s|(2s? - 1) + 8 s 2 c 2 c 3 cis 3 + 4c? (2s? - l)}s 8 
+S?{(8c?S? — 8c?S? + 8s? + 2c 2 — 8s?)s§ 

+8cic 2 s 2 s 3 (2c? - l)c 3 - 2c? (2s? - l) 2 }s 4 , (C.7i) 

0220 = ~Si s 3 s12 - S?{(2(2s? - 1 ))( c l + ^)4 - ( 2 ( 2s 2 ~ l))s§ + 8s?ClC 3 S 2 s|c 2 }s 8 
+{(—4c?s? + 4c?s? — 24c 2 s? — cf + 24c?s? — 4 s? — 2 cf + 4s? — l)s? 

+(24 44 - 24 44 + 8s? + 2 c? - 8s? + 2 )s| 

+8s 2 c 2 c 3 ci(2c? - l)(c?c§ - 4 + d)s 3 - (2s 2 2 - l) 2 }s 4 , (C.7j) 
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0211 = -2.sfc§s§s 12 + «i{—4(2c| - l)(c 2 + 1)4 + 16cis 2 c 2 s 3 c| 

+4(2c| - 1 )(c 2 + l)c§ - 8 cis 2 c 2 s 3 c 3 }s 8 

+{{&c \4 — 8 c \4 + 48 c 3 C2 + 2 c\ — 48c 2 c 2 + 8C2 + 4 c\ — 8c| + 2)c 3 
-16cic 2 s 2 s 3 (2c 2 - l)(c 2 + l)c| 

+(—8c^C2 + 8 c \4 — 48 c 2 cf — 2 c\ + 48 c 2 c| — 8C2 — 4 c 2 + 8 c| — 2)4 
+ 8 cic 2 s 2 s 3 ( 2 c| - l)(c 2 + 1 )c 3 - 8c|c 2 s|}.s 4 , (C. 7 k) 

0202 = -Si c 3 s12 + sf{—2(2 s\ - 1 )(c 2 + 1)4 - 8cis 2 c 2 s 3 c| + 2c§(2s| - l)}s 8 
+{(—4c 4 S2 + 4c 4 S2 — 24 c ^2 — c\ + 24 c ^2 — 4 s\ — 2 c\ + A.s^ — 1)03 
+8cic 2 s 2 s 3 (2c 2 - 1) (c 2 + l)c| + ( 24 c 2 ^ - 24 c ^2 + 8^ + 2c 2 - 8^2 + 2)c§ 
- 8 cic 2 s 2 s 3 ( 2 c| - l)c 3 - (2 s\ - l) 2 }s 4 . (C. 71 ) 


6423 = -s?s§, 6432 = -s\c\, bill = -Cl, (C.8a) 

^243 = Si{(- 4 c 2 S 2 C 2 - 4 s 2 c| + 1)4 ~ 4cic 2 s 2 s 3 (24 - l)c 3 - (2 si - l) 2 }s 4 , (C.8b) 

^234 = ■Si{(+4c 2 s 2 C2 + 4s\4 ~ 44 + 4cic 2 s 2 s 3 (2c2 - l)c 3 - 4 c 2 S2 c|}s 4 , (C.8c) 

6231 = {(4cjsi — 4 c 4 S 2 + 24 c\s\ — 24c 2 s 2 + 4 sf + 4c 2 — 4 s 2 )c|s| 

-4cic 2 s 2 S3(2c| - l)(c 2 + l)c 3 (l - 2 4 ) - 4( 2 4 - l) 2 }s 4 , (C.8d) 

6222 = {(—8c 4 C2 + 8 c \4 ~ 48 cfc 2 + 48 c 2 4 ~ 8cf — 8c 2 + 8 4)44 

+ 8 cic 2 s 2 s 3 (24 - l)(c 2 + l)c 3 (l - 2 4) 

+8c|(c 2 — l)(c 2 + l)(c 4 — c 2 + l)}s 4 , (C.8e) 

^213 = {( 4 cf 4 — 4 cfs| + 24 c 2 s| — 24c 2 s| + 4 + 4 c 2 — 4 s 2 )c|s| 

-4cic 2 s 2 S3(2c| - l)(c 2 + l)c 3 (l - 2 4 ) - c\(2sl - l) 2 }s 4 . (C.8f) 

^210 = 4{4444 + 4 c 2 s 2 - l)c§ + 4 cic 2 s 2 S3(2c| - l)c 3 - 4 c|c 2 s|}s 4 , (C.8g) 

b 2 01 = s 2 {( 4 c 2 44 ~ 444 + l)cl - 4cic 2 s 2 S3(2c 2 - l)c 3 - ( 2 c| - l) 2 }s 4 , (C.8h) 

^033 = -c?s 8 , 6012 = -S?c|s 8 , 6021 = -s 2 s|s 8 . (C.8i) 


C 442 = «i{(2ci - l)(c 2 + 1)4 - 4cic 2 s 2 S3c| + (—2c| + l)c|}s 4 - s 4 c 3 , (C.9a) 

C 433 = s?{8cic 2 s 2 s 3 c| + (2(2c 2 - l))(c 2 + 1)44 ~ 4s 2 c 2 c 3 cis 3 } s 4 - 24s\s\, (C.9b) 

C 424 = s?{-4cic 2 s 2 s 3 c 3 - (2 4 ~ 1)(2 4 + l )44 + 4s 2 c 2 c 3 cis 3 + c 2 (2c| - l)}s 4 

-Sl«3> ( C ' 9c ) 

C 421 = s 2 {-2c 2 (2c! - 1)4 + 4s 2 c 2 c 3 cis 3 }s a - 2.s 2 c 2 c§, (C.9d) 

C 412 = s?{-4s 2 c 2 c 3 cis 3 + 2c 2 s|(2s 2 - l)}s 4 - 2s 2 c 2 s|, (C.9e) 

C 400 = -sf4( 2 4 - !) s4 - c i> (C.9f) 

C244 = -sl4i 2 4 - i) 5 ® - 4 ( 2 4 - 1 ) 2 S 4 , (C.9g) 
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C 232 — 4{ — 4s 2 C 2 C 3 CiS 3 + 2cf(2s2 — l)s|}s 8 + 4{(^44 — 8 44 + 84 + 2 cf — 84)4 

-8cic 2 s 2 S3(2cl - l)c 3 - 2 c\{ 2 sl - l) 2 }s 4 , (C. 9 h) 

C 223 = 4 { 2 c i (24 — 1)4 + 4 s 2 c 2 C 3 C 1 S 3 } s 8 + sj{(— 844 + 8 44 — 8 s 4 — 2c 2 + 84)4 

+8cic 2 s 2 S3(2cl - l)c 3 + 8 s|(« 2 - 1)(S2 + l)}s 4 , (C. 9 i) 

C 220 = 4{( 2 4 - 4(4 + 1)4 - 4cis 2 c 2 s 3 c| 

— (2c 2 — l)( 2 c 2 + 1)4 + 4cis 2 c 2 s 3 c 3 + 4( 2 4 — l)}s 8 

+{(—4c 4 4 + 4c 4 c| — 24 c 2 4 — c\ + 2444 — 4cf — 2c 2 + 4c| — 1)4 

+8cic 2 s 2 s 3 (2c| - i)(c 2 +1)4 + 24(444 - 444 + i2 4 + 4 - 124 + 1)4 
-8s 2 c 2 cfs 3 (24 - l)c 3 - 4( 2 4 - l) 2 } s4 > (C. 9 j) 

c 2 n = 4{8cic 2 s 2 S34 - 2(24 - 1)(4 + 1)44 “ 4 s 2 c 2 c 3 cis 3 }s 8 

+{—16cic 2 s 2 s 3 (24 — l)( c i + 1)4 

+(—844 + 844 — 4844 — 24 + 4844 — 84 — 4 cf + 8s 2 — 2)44 
+8cic 2 s 2 s 3 (24 - 1)(4 + l) c 3 - 8444}s 4 , (C.9k) 

C 202 = 4{-( 2 4 -1)(4 +1)4 - 4cic 2 s 2 s 3 4 + ( 2 4 - i)c §}« 8 

+{(—444 + 4 44 — 244 4 — 4 + 2444 — 4 .s 2 — 24 + 4s 2 —1)4 

+8cic 2 s 2 s 3 (24 - 1)(4 + 1)4 + (2444 - 24 c\s\ + 8 4 + 24 - 8 s| + 2)4 

-8cic 2 s 2 s 3 (24 - l)c 3 - (24 - l) 2 }s 4 - (C.91) 
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